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I.  INTRODUCTION 

Prior  to  the  grant  period  5/15/82  -  5/14/83,  the  efforts  of  our  group 
concentrated  on  charge-exchange  schemes  to  obtain  lasing  in  the  X-ray  region. 
During  the  period  5/15/82  -  5/14/83,  the  efforts  of  the  group  were  directed 
towards  alternative  schemes,  such  as  the  free-electron  laser,  to  obtain  X-ray 
lasing.  Possible  uses  of  X-ray  lasers  in  novel  areas,  such  as  in  nuclear 
phenomena,  were  also  investigated,  Which  led  to  rigorous  analyses  of  some 
earlier  work  concerning  nuclear  beta  decay,  as  well  as  careful  considerations  of 
intense  field  behavior  of  quantum  systems. 

The  following  report  consists  of  two  parts.  The  first  part  summarizes  the 
initial  work  concerning  the  possibility  of  scaling  the  free-electron  laser  to 
X-ray  wavelengths.  The  second  part  summarizes  the  work  dealing  with  the 
possible  uses  of  X-ray  lasers  in  nuclear  phenomena,  and  with  correct 
descriptions  of  quantum  systems  in  intense  fields.  The  papers  which  detail 
these  results  either  in  published  or  in  preprint  form,  are  reproduced  as 
appendices . 

II .  SCALING  OF  FRLF-ELECTKON  LASER  TO  X-RAY  WAVELENGTHS 

A  difficulty  in  extending  the  operation  of  lasers  based  on  atomic 
transitions  down  to  X-ray  wavelengths  is  the  short  lifetime  of  excited  states, 
due  to  radiative  and  Auger  transitions.  This  difficulty  is  not  present  in  the 
case  of  free-electron  lasers  (FEL),  since  the  energy  source  for  the  lasing  is 
stored  stably  (for  hours,  in  the  case  of  storage  rings)  in  the  form  of  electron 
kinetic  energy.  Moreover,  the  FEL  is  intrinsically  a  swept-gain  device, 
normally  operating  with  overlapping  picosecond  pulses  of  electrons  and  laser 
radiation.  Thus,  a  gain  length  of  many  meters  is  possible  in  a  mostly  empty 
cavity,  since  the  gain  pulse  (electrons)  moves  at  nearly  the  same  speed  as  the 
laser  pulse.  The  FEL  wavelength  is  continuously  tunable  by  changing  the 
electron  energy,  and  so  can  be  matched  to  available  Bragg  crystal  resonators. 


Two  possible  ways  to  achieve  short-wavelength  operation  of  the  FEL  have 
been  considered.  One  is  to  use  high-energy  (GeV)  electrons  in  conjunction  with 
a  long  conventional  wiggler  (or  a  GW  microwave  field  in  a  waveguide).  The 
second  is  to  use  low-energy  (10  MeV)  electrons  passing  throu^i  a 
counter-propagating  high-power  infrared  laser  pulse.  The  first  regime  can  give 

t 

very  large  gain  (operation  without  a  resonator  is  a  possibility),  but  requires  a 
wiggler  tens  to  hundreds  of  meters  ion  length,  unless  a  good  quality  Bragg 
resonator  can  be  developed.  Electron  beam  quality  is  the  most  important 
limitation.  Storage-ring  operation  may  be  feasible  down  to  about  60  A.  Shorter 
wavelengths  (down  to  a  few  A)  are  possible  in  principle  using  a  linear 
accelerator,  but  suitable  accelerators  have  not  been  developed. 

The  second  regime  using  a  high-power  laser  pulse  can  give  moderate  (40%) 
gain.  In  this  regime  the  quantum  recoil  can  be  as  large  or  larger  than  the 
homogeneous  bandwidth,  so  that  quantum  effects  are  definitely  important.  The 
quantum  effects,  though  physically  interesting,  are  in  general  adverse,  reducing 
the  gain  and  limiting  the  saturation  power  to  about  one  photon  per  electron. 

The  correct  gain  can  be  calculated  using  semiclassical  theory  (quantizing  the 
electrons,  but  not  the  field).  Another  inportant  feature  of  this  type  of  FEL  is 
the  slowing  of  the  electrons  as  they  pass  through  the  high-intensity  beam  waist 
of  the  counterpropagating  Gaussian  beam.  This  slowing  is  due  to  transfer  of 
electron  kinetic  energy  into  transverse  motion  (mass-shift  effect)  and  spreads 
the  gain  over  a  bandwidth  dependent  on  the  pump  laser  power  (but  not  on  the  size 
of  the  beam  waist).  In  general,  electrons  in  an  appropriate  energy  range  are 
resonant  at  one  point  entering  and  at  one  point  leaving  the  focal  region,  and 
interference  between  radiation  emitted  at  these  two  points  gives  a  rapidly 
oscillating  homogeneous  spectral  lineshape.  The  highest  peak  of  the  lineshape 
curve  is  at  the  low-frequency  end  of  the  curve,  and  corresponds  to  electrons 
which  are  resonant  at  the  beam  waist.  The  frequency  of  this  peak  depends  on  the 
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power  of  the  pump  laser.  A  difficulty  with  this  regime  is  to  devise  a  cavity 
for  injecting  and  recirculating  the  terawatt  pump  pulse  and  to  keep  the 
electrons  resonant  for  many  passes  if  the  pump  pulse  is  decaying. 

The  work  sunsnarized  above  is  described  in  the  paper  by  J.  Gea-Banachloche , 

f 

G.  T.  Moore  and  M.  O.  Scully,  "Prospects  for  an  X-Ray  Free-Electron  Laser",  in 

,  i 

Free-Electron  Generators  of  Coherent  Radiation,  edited  by  C.  A.  Brau,  S.  F. 
Jacobs  and  M.  O.  Scully,  Proc.  SPIE  453,  393  (1984).  It  is  reproduced  in 
Appendix  I. 

III.  NUCLEAR  SPECTROSCOPY  WITH  X-RAY  LASERS 

One  can  reasonably  expect  that  X-ray  lasers  When  available,  will  have  a 
strong  impact  on  nuclear  spectroscopy,  similar  to  the  inpact  that  optical  lasers 
have  had  on  atomic  spectroscopy.  It  is  particularly  interesting  to  anticipate 
the  possibility  of  modulating  nuclear  beta  decay  or  orbital  electron  capture 
rates  by  means  of  an  X-ray  laser.  This  possibility  was  analyzed  for  a  simpl2 
case  of  a  parent  nucleus  (Z,N)  decaying  into  a  daughter  nucleus  (Z  ±  1,  N  t  1), 
where  the  parent  nucleus  has  just  two  states,  a  (lower)  and  b  (upper),  and  the 
daughter  nucleus  has  just  one  state  c.  Both  a  and  b  are  assumed  to  be  unstable 
with  respect  to  weak  decays.  The  X-ray  laser  frequency  is  assumed  to  be  tunable 
to  the  level  separations  of  a  and  b.  The  situations  similar  to  the  highly 
forbidden  decay  of  **|l  to  a  low  lying  excited  state  of  ^Xe.  physical 
picture  is  exactly  analogous  to  Raman  scattering.  The  results  of  the 
calculation  show  that  the  modulation  of  the  decay  rate  is  extremely  sensitive  to 
the  tuning  of  the  laser  frequency  to  the  a  **■  b  transition  frequency,  which  has 
natural  quantum  limits  arising  from  the  spontaneous  decay  of  the  states,  as  well 
as  from  internal  conversion  processes.  If  one  can  achieve  tuning  on  the  order 
of  Aw/w  *v  10'13,  powers  on  the  order  of  1014  W/cm2  are  needed  in  order  to  see 

observable  effects  in  the  decay  of  129 1 .  This  calculation  is  described  in  the 

53 

paper  by  W.  Becker,  R.  R.  Schlicher,  M.  0.  Scully,  M.  S.  Zubairy,  and  M. 


Goldhaber,  "Nuclear  Spectroscopy  with  X-ray  Lasers",  Ehys.  Lett.  13 1B(  1983)  16, 
which  is  reproduced  in  Appendix  II. 

One  interesting  result  concerns  scattering  processes  or  decay  of  charged  or 

neutral  pointlike  particles  in  a  laser  field  idealized  by  a  plane  wave  of 

y 

arbitrary  polarization.  The  wave  functions  of  the  charged  particles  are  then 

Volkov  solutions.  In  the  paper  by  W.  Becker,  G.  T.Moore,  R.  R.  Schlicher  and  M. 

O.  Scully,  “A  note  on  total  cross  sections  and  decay  rates  in  the  presence  of  a 

laser  field",  Rhys.  Lett.  94A  (1983)  131,  which  is  reproduced  in  Appendix  III, 

it  is  shown  that  in  the  quasiclassical  limit  (ft  -►  0),  the  differential  cross 

section  or  decay  rate  can  be  represented  as  an  integral  with  respect  to  time 

over  an  instantaneous  cross  section  of  decay  rate  Which  depends  on  the  value  of 

the  laser  field  at  a  given  time.  Since  this  instantaneous  rate  is,  in 

principle,  a  measurable  quantity,  it  must  depend  on  the  canonical  momentum  of 

the  ith  particle  and  on  the  laser  field  £(t)  only  via  the  mechanical  momentum 

*^(t)  =*  -  eiA(t).  If  the  differential  rate  is  integrated  over  the  final 

momenta  in  order  to  obtain  the  total  cross  section  of  the  total  decay  rate,  the 

dependence  of  the  final  rate  on  the  laser  field  is  completely  eliminated  by 

changing  the  integration  variable  from  ^  to  -  e.A(t).  Consequently,  the 

total  decay  rate  of  a  neutral  particle  is  unaffected  by  a  laser  field,  and  a 

charged  particle,  if  stable  in  vacuum,  remains  stable  in  a  laser  field.  The 

quasiclassical  approximation  breaks  down  if  the  strength  of  the  external  field 

becomes  comparable  with  the  critical  field  strength  E  =m2c3/eh'vl016  Volt/cm, 

crit 

which  is  out  of  reach  of  present  lasers  or  those  anticipated  in  the  near  future. 

The  preceding  considerations  can  be  extended  to  nucleons  bound  in  a 
self-consistent  potential  when  the  nucleon-laser  field  interaction  is 
negligible.  Given  the  energy  of  a  typical  laser  quantum  of  1  eV  and  a  typical 
nuclear  level  difference  of  1  MeV,  this  is  an  excellent  approximation.  In  this 
case  it  can  be  shown  again  that,  in  the  quasiclassical  limit,  total  decay  rates 
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are  unaffected  by  laser  fields  Which  are  weaker  than  the  critical  field.  Recent 
claims  that  forbidden  nuclear  beta  decay  can  be  drastically  enhanced  by  an 
intense  radio  frequency  field  are  contested.  This  work  is  described  in  Appendix 
IV,  which  is  a  preprint  entitled  "A  no-go  theorem  concerning  the  enhancement  of 
nuclear  decays  by  intense  radiation  fields”,  by  W.  Becker,  R.  R.  Schlicher  and 
M.  O.  Scully  (submitted  for  publication  in  Phys.  Lett.). 

At  the  1983  NATO  Advanced  Study  Institute  in  Madrid,  Spain,  M.  O.  Scully 
gave  seminars  on  the  subject  of  allowed  nuclear  beta  decay  in  intense  laser 
fields.  The  notes  for  these  lectures  were  prepared  by  R.  R.  Schlicher,  W. 

Becker  and  M.  O.  Scully  and  are  reproduced  in  Appendix  IV.  The  notes  cover 
topics  such  as  total  decay  rates  and  differential  electron  energy  distributions 
of  nuclear  beta  decay  in  the  presence  of  a  laser  field.  The  emphasis  is  on  the 
classical  versus  quantum  mechanical  description,  questions  of  gauge  invariance, 
and  the  fact  that  the  classical  description  is  sufficient  for  a  much  wider  range 
of  parameters  than  one  might  initially  assume. 
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Abstract 


He  present  an  overview  of  the  regimes  in  which  operation  of  an  x-ray  free-electron  laser 
(PEL)  may  be  feasible,  including  discussion  of  static  and  electromagnetic  wigglers,  quantum 
recoil,  high-gain  operation,  mass-shift  broadening,  and  electron  beam  quality. 

introduction 


The  free-electron  laser  is  a  swept-gain  laser  which  offers  unique  possibilities  for 
scaling  to  short-wavelength  operation.  Atomic  charge-exchange  swept-gain  devices1’2  which 
have  been  proposed  for  operation  at  x-ray  wavelengths  are  limited  by  rapid  spontaneous  decay 
of  the  upper  lasing  level.  By  contrast,  energy  input  to  the  PEL  is  in  the  stable  form  of 
electron  kinetic  energy,  and  the  electrons  themselves  (not  just  the  gain)  propagate  at 
nearly  the  speed  of  light.  NO  atomic  nuclei  or  bound  electrons  are  present  to  complicate 
the  physics  (except  in  the  resonator  materials)  . 

He  shall  consider  here  possible  regimes  for  an  x-ray  FEL  lasing  at  the  fundamental 
frequency  as  ■  2cys2k„  or  ■  4CYs2ki#  where  ys  is  the  longitudinal  Lorentz  factor  of  the 
electron  beam  and  kq  (k j)  is  the  wave  vector  of  the  static  (electromagnetic)  wiggler. 
Generation  of  higher  harmonics3 *4  of  u>s  is  another  important  mechanism  for  producing  short 
wavelengths,  but  will  not  be  treated  in  this  paper. 

Because  available  resonators  for  x-rays  are  of  poor  quality,  it  is  necessary  to  have 
large  gain  per  pass  to  obtain  lasing.  If  the  gain  is  sufficiently  large,  one  can  dispense 
with  the  resonator  altogether  and  lase  by  amplified  spontaneous  emission  (ASE) .  It  is  not 
known  how  much  gain  is  required  for  an  FEL  to  lase  by  ASE.  This  mode  of  operation  would 
have  advantages,  since  the  repetition  rate  of  electron  pulses  (and  wiggler  pulses  in  the 
case  of  a  counterpropagating  electromagnetic  wiggler  field)  could  be  arbitrarily  low.  There 
would  also  be  no  problem  of  resonator  alignment.  The  coherence  time  would  be  limited  to  the 
slippage  time  L/2yszc,  where  L  is  the  wiggler  length.  In  other  words,  the  coherence  in 
terms  of  x-ray  periods  could  not  exceed  the  number  of  wiggler  periods.  Probably  one  could 
not  do  much  better  than  this  using  a  low-Q  resonator. 

He  shall  consider  two  types  of  wiggler,  the  uniform  wiggler  and  the  Gaussian  beam.  The 
uniform  wiggler  could  be  a  conventional  permanent-magnet  or  bif ilar-helical  type.  It  could 
also  be  a  counterpropagating  electromagnetic  wave  contained  in  a  wave  guide.  The  latter 
might  be  less  expensive,  particularly  if  one  needs  a  very  long  wiggler.  In  either  case  the 
wavelength  x  or  x  «  2x  is  on  the  order  of  centimeters.  To  obtain  x-rays  from  such  a 
wiggler,  the^electron  energy  must  be  in  the  GeV  range. 

The  Gaussian-beam  wiggler  uses  a  counterpropagating  focused  Gaussian  beam  from  a 
high-power  infrared  laser  as  the  wiggler  field.  Because  Xj  is  small,  the  electron  energy 
need  only  be  on  the  order  of  tens  of  MeV.  The  interaction  length  is  limited  both  by  the 
Rayleigh  range  2q  and  the  infrared  pulse  length.  Moreover,  slowing  of  the  electrons  in  the 
vicinity  of  the  beam  waist  is  the  dominant  source  of  homogeneous  line  broadening  for  a 
powerful  (terawatt)  wiggler  field.  Quantum  recoil  of  the  electrons  is  also  typically 
important  for  this  type  of  FEL. 


Low-Gain  Uniform  Higgler 

Let  us  now  consider  the  uniform  wiggler  in  detail.  The  quantity  of  most  interest  is  the 
small-signal  gain.  We  first  evaluate  the  gain  per  pass  G  assuming  G  <<  1.  Then  we 
generalize  to  the  large-gain  regime  G  >>  1.  For  the  sake  of  definiteness,  we  assume  a 
static  wiggler.  We  neglect  pulse  effects,  assume  the  x-ray  laser  beam  to  be  monochromatic, 
and  account  for  the  detuning  of  electron  energies  E  from  resonance  by  the  usual  detuning 
parameter  defined  by 

E  •  Mc2vs(l  ♦  u/21  )  .  (1) 


Here  m  •  m  j'  is  the  electron  mass  multiplied  by  a  mass-shift  correction  given  by 
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I  «  1  ♦  (eB/mck^) 2  ,  (2) 

where  B  is  the  RMS  magnetic  field.  In  practice  we  are  interested  in  cases  where  i  is  at 
most  a  few  times  unity. 

If  we  postpone  for  now  the  question  of  emittance,  there  are  three  detuning  widths 
governing  the  expression  for  the  gain.  First,  there  is  the  homogeneous  transit-time 
broadening  with  bandwidth  2n/L.  Second,  there  is  the  inhomogeneous  broadening  U 
characterizing  the  width  of  the  normalized  electron  energy  distribution  f(i-).  Third,  there 
is  the  quantum  recoil  2q  undergone  by  electrons  when  they  emit  an  x-ray  photon. 

The  scaling  of  U  is  different  for  linear  and  circular  accelerators  (or  storage  rings) . 

For  linear  accelerators  the  spread  4E  is  determined  mainly  in  the  initial  part  of  the 
accelerator.  With  care  one  can  accelerate  the  electrons  to  high  energy  without  increasing 
BE.  From  Eq.  (1)  we  see  that 

U  «  -2-k-3--  4E  -  E2--  6E  A  **  1  ',s  .  (3) 

Mc2y#  Me2  5  ** 

The  value  of  6E  depends  on  the  particular  accelerator,  but  probably  the  minimum  value  one 
could  obtain  is5 

SE  =  el/4rre0c  =»  mc2(I/l7,000  A)  .  .  (4) 

where  1  is  the  electron  current.  This  limit  comes  from  the  Coulomb  repulsion  of  the 
electrons.  The  energy  spread  in  the  Stanford  superconducting  LltlAC  is  about  100  times 
larger . 

In  a  storage  ring  the  dominant  energy  spread5  is  due  to  noise  associated  with  synchrotron 
radiation  and  gets  larger  at  higher  energies.  One  has  approximately 

SE  •  (h/2mco)!i(E2/mc2)  ,  (5) 

where  o  is  the  bending  radius  of  the  ring.  Combining  Eqs.  (5)  and  (1)  yields 

U  »  ZnCha/mco)*5  A#"^  A^""1  ,  (6) 

The  fact  that  Eq.  (6)  scales  as  suggests  that  a  storage-ring  x-ray  FEL  is  impracticable 

if  \s  is  too  small,  since  the  maximum  L  for  homogeneously  broadened  operation  will  scale  as 
As**.  Our  estimates  indicate  a  minimum  A$  of  about  60  A.  The  storage  ring  seems  an  excel¬ 
lent  option  for  A,  of  100  A  or  greater.7  With  a  linear  accelerator  it  may  be  possible  to 
scale  A$  down  to  5  A  or  less. 

The  quantum  recoil  when  expressed  in  detuning  units  becomes 

2q  «  (8  ir2  2^h/Mc)  A  A  .  ,7) 

s  q 

In  the  small-gain,  small-signal  regime  the  gain  may  be  written  as  G  «  C.p,  where 

c  .  -,,c3ib-*-..  A  *  .  (6) 
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We  have  used  the  ■'~«ole  upshift  condition  to  eliminate  explicit  reference  to  -.in  Eq.  (8). 
I  is  the  current  wjt  in  the  laser  mode  area  r.  The  coefficient  r(  is  one  for  a  helical 
wiggler  and  a  well  .nown 3  difference  of  Bessel  functions  for  a  linear  wiggler.  Equation  (9) 


expresses  r  as  the  difference  between  the  forward  and  reverse  quantum  processes.  This 
difference  is  replaced  by  a  derivative  with  respect  to  u  in  classical  treatments  of  the  FEL. 
This  replacement  is  sometimes,  but  not  always,  justified  for  the  x-ray  regime. 

We  may  distinguish  six  regimes  depending  on  the  relative  sizes  of  the  homogeneous, 
inhomogeneous,  and  quantum  widths,  as  shown  in  Table  1. 

As  an  example,  let  ufs  consider  a  helical  wiggler  with  the  parameters  of  the  Stanford 
wiggler,  except  that  it  is  longer,  and  use  this  wiggler  to  generate  5  A  radiation.  The 
system  parameters  we  use  are  shown  in  Table  2. 


2q  <  U  <  2x/L, 

R 

2q  <  2x/L  <  0, 

n 

0  <  2q  <  2x/L, 

2  x/L  <  2q  <  0, 

0  <  2x/L  <  2q, 

2 x/L  <  0  <  2q, 

■ 

Table  1 


Quantity 

Symbol 

Value 

Wiggler  wavelength 

3.2  cm 

Magnetic  field 

B4 

.24  T 

Mass  shift 

a 

1.512 

Laser  wavelength 

i 

5  A 

Electron  energy 

ES 

3.554  GeV 

Current 

I 

10  A 

Electron  beam  power 

El/e 

3.564  x  10!,3W 

Laser  mode  area 

r 

.  1 6  mm- 

Table  2 

With  these  parameters  we  calculate  a  homogeneously  broadened  gain  of  unity  at  L  •  61  m.  If 
we  take  the  inhomogeneous  broadening  to  be  given  by  Eqs.  (3)  and  (4),  and  we  use  Eq.  (7)  for 
the  quantum  recoil,  we  calculate  2x/u  *  1.9  x  105m  and  x/q  =  2.3  x  10'm.  These  correspond 
to  the  wiggler  lengths  at  which  the  homogeneous  broadening  would  equal  the  inhomogeneous 
broadening  or  the  quantum  recoil  respectively.  We  conclude  that  5E  could  be  more  than  three 
orders  of  magnitude  larger  than  the  value  given  by  Eq.  (4)  before  inhomogeneous  broadening 
begins  to  affect  the  gain. 

So  far  we  have  merely  assumed  that  \  *  3.2  cm  without  justifying  whether  this  is  close 

to  an  optimal  value.  We  can  get  a  more'‘global  picture  of  the  situation  on  a  log-log  plot  of 

L  versus  ,\q,  as  shown  in  Fig.  1.  Here  the  point  of  unit  gain  at  .\-  =  3.2  cm,  L  *  81  m  is 

marked  by  an  X.  The  two  lines  with  slope  3/2  indicate  the  transitions  where  -/L  =  q  and 

2x/L  =  U.  If  we  hold  B  and  Z  constant,  and  neglect  minor  effects  due  to  mass  shift,  we  can 
generate  a  curve  along  which  G  =  1  by  using  Eq.  (8)  and  Table  1.  This  gives  the  three- 
segment  solid  curve  in  Fig.  1.  Alternatively,  if  we  assume  that  the  mode  area  r  is  governed 
by  diffraction  of  the  x-rays,  then  z  «  1SL  and  we  generate  the  dotted  curve  in  Fig.  1.  We 
can  generate  curves  along  which  the  gain  has  other  values  by  parallel  transport  of  these 

curves  along  the  transition  lines.  These  curves  of  constant  gain  are  invalid  if  1  becomes 

too  large.  This  is  because  the  mass-shift  corrections  then  become  substantial  and'because 
the  wiggle  amplitude  may  exceed  the  assumed  mode  size. 

We  can  infer  from  Fig.  1  that  it  is  best  to  choose  conditions  close  to  the  transition 
lines,  since  L  is  minimized  there.  The  choice  >q  =  3.2  cm  gives  us  some  latitude  if  :E  is 

larger  or  if  we  want  to  increase  L  further  to  get  more  gain.  We  also  see  that  it  is  very 

disadvantageous  to  decrease  '  to  operate  in  the  inhomogeneously  broadened  regime.  The  main 
point  in  trying  to  use  a  terawatt  infrared  laser  as  a  wiggler  is  simply  that  the  magnetic 
field  one  can  generate  is  much  larger  than  in  a  conventional  wiggler. 

High-Gain  Uniform  Wiggler 

The  classical  cold-beam  small-signal  equations  for  the  FEL  have  been  published 
previously.5  The  extension  of  these  equations  to  include  inhomogeneous  broadening  is 
straightforward.  The  extension  to  include  quantum  recoil  is  too  lengthy  to  derive  here,  but 
can  be  done  by  analyzing  the  coupled  Maxwell  and  Xlein-Gordon  equations.  The  result, 
itself,  is  very  simple.  In  the  cold-beam  limit  the  equations  take  the  form 

elB'\  "■  -iuo- 
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Figure  1.  Curves  of  constant  gain  on  a 
log-log  plot  of  wiggler  length  L  versus 
wiggler  wavelength  1  provide  a  convenient 
global  picture  of  regions  where  operation  of 
a  short-wavelength  FEt.  may  be  possible.  In 
this  figure  we  fix  the  laser  wavelength,  the 
wiggler  magnetic  field,  and  the  current  at 
the  values  given  in  Table  2.  The  point 
marked  X  corresponds  to  the  unit-gain 
example  discussed  in  the  text.  The  energy 
spread  is  held  constant  at  300eV, 
independent  of  the  electron  energy.  The 
quantum  recoil  in  this  example  exceeds  the 
inhomogeneous  broadening,  and  the  two 
parallel  lines  of  slope  3/2  indicate  the 
transitions  where  these  respective  widths 
equal  the  homogeneous  broadening.  Using  Eq. 
(8)  and  Table  1,  and  assuming  the  laser  mode 
area  \  is  constant,  we  construct  the 
three-segment  solid  curve  along  which  G  »  1. 
The  slopes  of  the  three  segments  (from  right 
to  left)  are  1/6,  -1/2  and  -5/2.  By 
parallel  transport  of  this  curve  along  the 
transition  lines,  we  can  generate  curves 
along  which  G  has  other  values.  If  G  is 
small,  then  G*A  4  at  the  transitions.  The 
value  of  G  on  a^given  curve  scales  as  B2I/i. 
We  see  that  it  is  optimal  (in  terms  of 
minimizing  L)  to  operate  near  the  transition 
lines.  Also,  it  is  disadvantageous  to 
choose  A  too  small,  unless  one  can 
compensate  by  increasing  the  wiggler  field. 

\  «  asl,  as  is  appropriate  if  the  mode  area 
this  case  the  segments  of  the  curves  of  cons 
gain  along  the  transition  lines  scales  as  \q 
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The  dotted  curve  is  generated  by  taking 
S  limited  by  diffraction  of  the  x-ray  oeam.  Ir 
ant  gain  have  slopes  1/4,  -1  and  and  t.ne 
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Here  Es,  K2,  and  Ki  are  the  laser  field,  the  density-bunching  amplitude,  and  the 
velocity-bunching  amplitude.  The  last  term  in  Eq.  (13)  is  the  quantum  correction,  and  is 
the  velocity  detuning  from  resonance.  The  quantum  correction  arises  from  an  interference 
between  probability  amplitudes  with  and  without  the  photon  emission.  In  the  small-gam 
limit  Eqs.  (11)— (13)  can  be  used  to  derive  (8)— (10)  in  the  case  f(u)  »  5(u  -  -j).  However, 
in  general  the  solution  of  Eqs.  (11)— (13)  consists  of  three  exponential  modes  having  the 
spatial  dependence  exp  (iz),  where  6  obeys  the  cubic  dispersion  relation 


B[(3  *  iuo  ♦  U)2  ♦  q2]  -  iC 


In  Eq.  (14)  we  have  generalized  to  the  case  where  f(u)  is  not  a  5-function,  but  is  a 
Lorentzian  of  width  U  (HWHM)  centered  at  u  *  ^9.  At  most  one  of  the  roots  of  Eq.  (14)  has  a 
positive  real  part,  and  the  corresponding  mode  will  dominate  near  the  end  of  a  sufficiently 
long  wiggler.  At  the  entrance  of  the  wiggler  the  three  modes  contribute  about  equally  in 
order  to  satisfy  the  initial  conditions  K,(0)  «  K2(0)  =  0.  Figure  2  shows  the  growth  rate 
2  Re(d)  for  q  =  0,  U  =  0  as  a  function  of  yj.  We  see  that  maximum  gain  is  obtained  on 
resonance,  where 


1  I  1  •/, 

g  exp ( 3  2  C'5  L) 


For  the  parameters  of  Table  2  we  calculate  a  g3in  per  pass  of  1000  if  L  *  276  m.  This  might 
be  enough  to  lose  without  a  resonator.  A  conventional  wiggler  of  this  lenqth  would  oe  very 
expensive,  but  one  might  be  able  to  instead  use  a  microwave  confined  by  a  waveguide."  The 
microwave  power  required  would  be  about  1  GW,  and  the  minimum  microwave  pulse  length  needed 
would  be  2L  *  552  m,  so  that  the  microwave  pulse  energy  would  be  about  1.8  kd.  Sucn 
microwave  powers  are  obtainable  and  can  be  propagated  in  a  single  wave-guide  rode. 

Microwave  pulse  lengths  obtained  so  far  at  this  power  are  a  bit  short  of  what  we  would  like. 
Note  that  the  slippage  time  L/2rs‘c  is  only  1.44  x  10  '“sec. 


Figure  2.  The  gain  coefficient  2  Re(3) 
for  the  exponentially  growing  mode  is 
shown  as  a  function  of  velocity  detuning 
uo>  neglecting  quantum  recoil  and 
inhomogeneous  broadening. 


Figure  3.  The  gain  coefficient  2  Re(3) 
is  shown  as  a  function  of  inhomogeneous 
broadening  U  for  the  optimal  detuning 
m o  “  U/31*  Quantum  recoil  is  neglected. 


The  effect  of  a  Lorentzian  inhomogeneous  broadening  of  width  U  (again  neglecting  q)  is  to 
shift  the  optimal  detuning  to  uo  •  3  :U  and  to  reduce  the  gain,  as  shown  in  Fig.  3. 

In  Fig.  4  we  show  2  Re ( 8 )  for  the  case  where  the  inhomogeneous  broadening  U  is  set  equal 
to  zero,  but  the  quantum  recoil  q  is  taken  into  account.  If  C  J  q3 ,  then  the  gain  is  peaked 
near  i.0  «  q.  If  C  becomes  larger  than  q3»  then  2  Re(i)  approaches  the  classical  result 
given  in  Fig.  2,  but  continues  to  have  a  sharp  cut-off  on  the  negative,  as  well  as  positive, 
end  of  the  detuning  range.  The  quantum  recoil  reduces  the  maximum  of  2  Re(£)  below  its 
classical  value.  For  the  parameters  of  Table  2  this  is  a  very  slight  effect,  since 


2.7  x  10°. 


Figure  4.  The  gain  exponent  (in 
units  of  cv*)  vs.  the  detuning 
(in  units  of  q)  for  several 
cases  in  the  quantum-mechanical 
regime.  Solid  Line:  c/q3  ■  1. 
Dashed  line:  C/q3  »  .1. 

Dash-dot  line:  C/q*  =  .01. 


2  Re (0) 

C'/3 


We  expect  the  FEL  to  saturate  when  the  density  bunching  amplitude  K2(z)  becomes  of  the  order 
of  unity.  From  this  condition  we  can  derive  the  saturation  power 
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For  the  above  numerical  example  we  calculate  PS3t  =1.7  MW. 


Gauss ian-Beam  Wigqler 


Infrared  lasers  such  as  COi  or  tld-glass  can  generate  very  large  electromaonetic  fields, 
which  suggests  that  we  consider  high-power  pulses  from  sucn  lasers  as  wiggler  lields  for  an 
x-ray  FEL.*'  Since  is  short,  one  needs  relatively  low  electron  energies  (tens  of  MeV) . 
The  interaction  length  is  limited  to  half  the  pulse  length  of  the  infrared  (pump)  laser, 
since  the  electrons  and  pump  field  are  counterpropagating,  each  at  nearly  the  speed  of 
light.  Moreover,  the  interaction  length  is  limited  to  a  value  on  the  order  of  the  Rayleigh 
range  Z  by  diffractive  spreading  of  the  pump  field.  In  the  present  analysis  we  take  the 
pulse  length  infinite  for  simplicity  and  only  take  into  account  the  diffractive  spreading. 
Then  the  on-axis  pump  field  is  of  the  form 
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The  wiggler  field  is  therefore  tapered .  The  electrons  will  be  slowed  up  in  the  vicinity  of 
the  focus,  where  more  of  their  energy  gets  transferred  into  transverse  oscillation.  This 
results  in  greater  hompgeneous  broadening  than  for  a  uniform  wiggler.  Electrons  within  the 
homogeneous  bandwidth  will  in  general  be  resonant  (move  at  the  same  speed  as  the 
ponderomotive  potential)  only  at  one  point  entering  and  one  point  leaving  the  focal  region. 
The  electron  slowing  is  a  consequence  of  the  variable  mass  shift.  Even  though  I  -  1  will 
normally  be  very  small  for  the  Gaussian-beam  wiggler,  the  cumulative  effect  over  the  many 
wiggler  wavelengths  within  the  interaction  region  can  be  large.  It  turns  out  that  the  size 
of  the  effect  is  characterized  by  a  parameter  e  depending  only  on  the  power  p,  in  the 
pump  f ield ,  1 


c  ■  (2e2/ne om2cs) . 


For  a  terawatt  beam  this  gives  c  ■  913.  Let  us  specify  that  the  energy-detuning  u  be  zero 
for  electrons  which  at  z  ■  ±«  travel  at  the  speed  c(ks  -  ki)/(ks  +  ki)  of  the  ponderomotive 
potential.  Then  the  homogeneous  bandwidth  for  large  e  is  approximately  0  <  uZ„  <  e.  The 
function  a(w)  in  Eq.  (10)  is  replaced  by 


1  -  iz/Z. 


exp[-iwz  ♦  it  arctan(:/;R) ] 


Also,  the  constant  C  in  Eq.  (8)  can  be  written  in  terms  of  c  a: 


2c3e0i:n>Ys  3’p 


The  two  points  of  stationary  phase  in  Eq.  (19)  lie  at  z  *  s  ZR(c/uZR  -  l)"3,  and  A  is 
approximately  a  sum  of  contributions  from  these  two  points.  The  interference  of  these  two 
contributions  makes  |  A <  u )  I  2  a  rapidly  oscillating  function  of  j  if  £  is  large.  A  graph  of 
this  function  for  s  «  92  is  shown  in  Fig.  5.  The  width  of  the  large  peak  on  the  right  side 
is  approximately  2eVi/ZR.  It  is  possible  to  express  A ( u )  exactly  in  terms  of  Whittaker 
functions.11  An  approximate  expression  in  terms  of  Airy  functions  which  is  valid  near 
u  ■  e/ZR  (that  is,  around  the  main  peak  on  the  right  side)  is 
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Figure  5.  The  mass-shift  broadened 
line-shape  r(u)  as  a  function  of 
energy  detuning.  THe  main  peak 
on  the  right  (centered  at  about 
u  *  c/ZR)  corresponds  to  a 
detuning  such  that  particles  are 
resonant  near  the  focus  of  the 
Gaussian  beam. 
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where  Ai(x)  is  the  Airy  function, 


As  j  becomes  less  than  r-/Zn,  the  function  |'(u)|‘-  becomes  more  and  more  rapidly 
oscillating.  In  cases  where  the  inhomogeneous  broadening  is  large  compared  to  the  period  of 
these  oscillations,  we  can  average  over  the  oscillations  to  obtain  (in  the  stationary-phase 
approximation) 
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Contributions  to  } A ( u ) | 2  near  u  ■  0  come  from  z  points  far  away  from  the  beam  waist.  In 

practice  such  contributions  are  excluded  by  the  finite  pump  pulse  length. 

Clearly  there  are  a  number  of  possible  operating  regimes,  depending  on  the  relative  sizes 
of  U  and  2q,  compared  to  2c’’/Zr  and  e/ZR.  The  most  favorable  case  is  one  in  which  all  the 
electron  energies  are  contained  within  the  main  peak  on  the  right  of  the  gain  curve  (see 

Fig.  5) ,  and  where  the  Compton  recoil  is  of  the  same  order  as  the  width  of  that  peak,  so 

that  the  absorption  process  is  well  separated  from  the  emission  process  (with  no  electrons 
contributing  to  the  former).  In  other  words,  we  require 

U  <  4^  *  2q  .  (23) 

-R 

where  for  an  electromagnetic  pump  pulse  of  wavelength  1^,  2q  is  given  by 
„  32*2h  ,  -V4  , 

2q  *  -  X  X.  (24) 

M  me  s  l 

(This  corresponds  to  Eq.  (7)  with  X.  ■  2x  ).  Under  these  conditions,  using  Eq.  (21),  we 
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max  2c2c0S*  (26) 

where  we  used  Eq.  (20)  and  the  fact  that  the  maximum  value  of  Ai 2  (x)  is  0.29. 

Once  the  wavelength  of  the  pump  is  chosen,  the  value  of  the  ratio  e%/ZR  is  limited  by  I 
the  condition  (23).  Recall  that  c  is  proportional  to  the  power  in  the  pump  pulse  and  ZR  is 
proportional  to  the  pulse  length.  Increasing  ZR  helps  to  increase  the  gain,  according  to 
Eq.  (26)  ,  but  it  also  reduces  the  width  of  the  gain  peak,  and  eventually  one  ends  up  in  the 
inhomogeneously  broadened  regime. 

A  serious  limitation  to  the  minimum  achievable  energy  spread  comes  from  Eq.  (4),  which 
gives  the  difference  in  the  kinetic  energy  of  the  electrons  at  the  edge  and  the  electrons  on 
the  axis  of  an  unneutralized  beam.  The  inhomogeneous  width  resulting  from  this  effect  is 
given  by 
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Note  that  I  here  is  not  necessarily  the  total  current,  but  only  the  current  through  the 
laser  mode  area  I.  It  seems  as  if  one  could  make  this  very  small  by  making  the  laser  mode 
area  very  small  (hence  decreasing  I  while  keeping  l/l  constant).  In  reality  one  is  limited 
here  too  by  the  fact  that  the  electrons  have  some  transverse  velocity  and  if  I  is  very 
small,  they  might  drift  out  of  the  interaction  volume  and  contribute  little  to  the  gain.  A 
measure  of  this  transverse  velocity  spread  is  provided  by  the  emittance  of  the  beam. 

»  m 

Without  dwelling  on  the  details,  we  present  in  Table  3  some  optimal  values  for  two 
possible  cases:  one  in  which  the  pump  pulse  is  a  1012  W  pulse  from  a  Ud : YAG  laser  at  1.1 
i>m,  and  one  in  which  it  is  a  10' 1  w  pulse  from  a  CO*  laser  at  10.6  urn.  In  both  cases  we 
assume  operation  at  x  ■  5.7  K.  We  have  chosen  this  wavelength  because  of  the  possibility 
at  this  wavelength  or  building  a  reasonably  good  resonator  using  Bragg  diffraction  from  Ge 
crystals.10  We  take  the  length  of  the  pump  pulse  to  be  roughly  equal  to  the  Rayleigh  length 
to  obtain  efficient  coupling. 

The  C02-laser-pumped  case  appears  to  be  more  favorable,  except  that  the  relative  energy 
spread  is  rather  small.  The  absolute  energy  spread,  however,  is  the  same  in  both  cases,  and 


Parameters 


Pump  wavelength  a^ 

Pump  power  Pi 

Rayleigh  length  ZR 

Total  pulse  energy  ri(2Zn/c) 

1.1  um 

10 12  watts 
(e  -  913)  . 

4.4  mm 

30  Joules 

10«6  um 

10u  watts 
(c  -  91.3) 

6  cm 

40  Joules 

Electron  Beam  Parameters 

Electron  beam  energy  E 

11  MeV 

34  MeV 

Normalized  emittance 

5.2  x  lO'S  m  rad 

same 

Peak  current  density  1/  I 

10  M A/cm  2 

1 . 6  MA/cm  2 

Peak  current  (total) 

9  kA 

4.7  kA 

Energy  spread  6E/E 

2  x  10‘4 

6  x  10"5 

Electron  beam  cross  section 

9  mm2 

29  mra2 

Laser  mode  area 

7xl0"l°  m2 

4.4xl0*9m 2 

Gain 

0.4 

0.4 

Table  3 

equals  2.4  keV  (equal  to  the  Compton  recoil  at  the  wavelength  considered).  The  normalized 
emittance  has  been  chosen  equal  to  that  of  the  Stanford  FEL  experiment.  A  smaller  emittance 
would  help  lower  the  requirements  on  the  total  current.  Note  in  this  respect  that  the  total 
peak  current  across  the  laser  mode  area  is  in  both  cases  only  about  70  A,  so  that  many  I 
electrons  are  being  "wasted".  As  mentioned  before*  this  is  because  one  needs  to  keep  the 
electron  beam  area  large  so  that  it  won't  spread  too  fast,  given  the  emittance  we  have 
assumed.  The  need  for  high  current  densities,  on  the  other  hand,  seems  unavoidable,  as  may 
be  seen  from  Eq.  (26) . 

The  situation  is,  as  one  might  expect,  more  favorable  for  longer  As  a  final  example, 

we  give  some  figures  for  operation  at  100  A,  using  again  a  terawatt  pulse  from  a  Nd-YAG 
laser.  In  this  case,  we  increase  ZR  over  the  value  given  in  the  previous  example  to  satisfy 
Eq.  (23).  With  the  value  Z_  *  1.74  cm  (corresponding  to  a  total  pump  pulse  energy  of  120  J) 
one  gets  a  gain  of  unity  (lOOt)  for  a  current  density  of  21  kA/cm2  and  a  peak  current  of  190 
A.  The  electron  energy  is  only  2.6  MeV  in  this  example,  and  the  maximum  allowable  energy 
spread  to  stay  in  the  homogeneously  broadened  regime  is  4E/E  »  5  x  10'S  («E  •  0.1  keV) . 

Discussion 

In  this  paper  we  have  had  space  only  to  describe  the  basic  physics  of  various  regimes  for 
possible  operation  of  an  x-ray  FEL.  We  have  seen  that  substantial  gain  can  be  achieved. 
However,  clearly  much  work  remains  in  areas  such  as  accelerator  design,  electron  beam 
transport,  x-ray  resonators,  and  optics  for  injection  and  recirculation  of  high-power 
electromagnetic  wiggler  pulses.  The  subject  is  relatively  undeveloped,  and  additional 
regimes  of  operation  and  mechanisms  for  gain  enhancement  may  remain  to  be  discovered. 

We  conclude  by  assessing  several  items  of  concern  to  the  designer  of  an  x-ray  FEL. 
Diffractive  spreading  of  the  x-ray  beam  is  a  minor  effect,  although  appreciable  for  a  very 

long  system.  The  value  of  £  which  we  used  for  the  long-wiggler  example  is  consistent  with 

the  effect  of  diffractive  spreading.  For  the  Gaussian-beam  wiggler  one  can  have  a  very 
small  £.  However,  if  the  electron  beam  area  is  much  larger  than  £,  it  could  be  hard 
experimentally  to  observe  the  lasing  against  a  strong  background  of  incoherent  emission  from 
all  the  electrons.  Because  of  the  high  loss  of  x-ray  resonators,  one  can  expect  the  laser 
output  itself  to  be  rather  incoherent. 

Absolute  tolerances  on  the  long  conventional  (or  microwave)  wiggler  are  about  the  same  as 
in  present-day  conventional  wigglers./  For  the  Gaussian-beam  wiggler  the  pump  beam  must  be 
coherent  over  the  interaction  length.  Effects  due  to  the  finite  pulse  length  of  the  wiggler 
pulse  may  also  be  important,  and  remain  to  be  calculated.  If  one  operates  with  £  >>  l,  then 

either  the  power  in  the  pump  pulse  must  be  kept  nearly  the  same  on  each  pass  or  the  electron 

energy  must  vary  from  pass  to  pass  in  order  to  maintain  the  gain  at  constant  x-ray 
wavelength . 

Transverse  variations  in  the  wiggler  field  are  a  serious  concern  for  all  the  regimes  we 
have  considered,  since  they  can  lead  to  de-phasing  of  the  electrons  with  respect  to  the 
ponderomotive  potential.  We  estimate  that  such  effects  can  be  adequately  controlled  by 
keeping  the  e-beam  area  small  enough  that  the  electrons  essentially  see  the  on-axis  wiggler 
field. 

Emittance  can  also  reduce  the  gain.  We  estimate  that  a  normalized  emittance  equal  to  or 
less  than  that  of  the  Stanford  superconducting  accelerator  is  sufficient  to  prevent  serious 
gain  depletion. 


t 


One  reason  why  it  is  interesting  to  study  the  physics  of  the  x-ray  FEL  is  that  features 
of  the  problem  which  are  of  minor  importance  at  long  wavelengths  here  become  dominant.  In 
particular,  quantum  and  other  fluctuations,  as  well  as  the  development  of  coherence  in  three 
-  -dimensions  through  ASE,  are  dominant  processes.  Future  theoretical  work,  together  with 
experiments  at  short  wavelengths,  should  lead  to  a  better  understanding  of  this  largely 
uncharted  territory. 
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The  effects  expected  from  irradiating  a  beta-unstable  nucleus  having  a  very  low-lying  excited  state  with  a  resonant 
X-ray  laser  are  investigated. 

Much  as  the  availability  of  optical  lasers 
brought  about  a  revolution  in  atomic  spectros¬ 
copy,  strong  short  wave  length  X-ray  or  gamma- 
lasers  would  induce  a  comparable  breakthrough 
in  nuclear  spectroscopy.  In  this  note  we  want  to 
determine  the  laser  specifications  which  would 
be  necessary  for  that  purpose.  As  an  example 
we  shall  treat  a  situation  as  depicted  in  (ig.  1.  A 
parent  nucleus  (Z,  N)  in  its  ground  state  a  is 
assumed  to  decay  via  beta*-decay  or  orbital 
electron  capture  to  the  state  c  of  the  daughter 
nucleus  (Z  s  I,  N  *  1)  which  can  be  the  ground 
state  or  an  excited  state.  In  the  presence  of  an 
incident  laser  field  which  is  closely  resonant 
with  the  excitation  energy  of  an  excited  level  b 
of  the  parent  nucleus  the  decay  can  proceed 
alternatively  via  absorption  of  one  laser  photon 
and  subsequent  beta-decay  from  the  level  b.  In 
the  latter  case  the  beta-decay  is  governed  by  the 
weak  interaction  matrix  element  V*  instead  of 
Vjc  which  specifies  the  ordinary  beta-decay.  By 
comparing  the  decays  with  and  without  the  laser 
field,  information  can  be  obtained  about  the 
matrix  element  V*  which  is  inaccessible  other¬ 
wise.  In  particular,  if  the  direct  decay  from  a  to 
c  is  forbidden  whereas  the  decay  from  b  to  c  is 
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Research  under  grant  No.  AFOSR-8I-OI28. 

:  Supported  by  the  US  Department  of  Energy  under  Con¬ 
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(Z.N)  (Zal.Ntl) 

Fig.  1.  The  nuclear  energy-level  configuration,  which  is 
envisioned  in  the  present  paper;  the  parent  nucleus  (Z.  N) 
with  ground  state  a  has  a  low-lying  excited  state  b  which  is 
almost  resonant  with  an  incident  X-ray  laser  of  frequency  to. 
Both  states  a  and  b  are  beta-unstable. 

allowed,  a  large  enhancement  of  the  total  decay 
can  be  achieved.  A  similar  situation  has  been 
referred  to  as  excited  state  beta-decay  in  astro¬ 
physics  [1].  However,  in  the  case  we  are  con¬ 
sidering,  the  decay  proceeds  virtually  via  the 
intermediate  state  b,  whereas  at  temperatures 
considered  in  astrophysics  the  intermediate  state 
is  thermally  populated.  Also,  a  similar  level 
scheme  has  been  suggested  in  order  to  populate 
the  upper  level  (level  b  in  our  notation)  of  a 
possible  nuclear  X-ray  laser  (2). 

Our  formal  approach  is  exactly  analogous  to 
the  treatment  of  Raman  scattering.  We  assume 
that  the  laser  is  almost  resonant  with  the  energy 
separation  from  a  to  b.  The  Schrodineer  equa¬ 
tion  then  yields  a  system  of  equations  for  the 
amplitudes  of  the  levels  a,  b  and  c.  In  the  in¬ 
teraction  representation  it  reads  [3] 
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id  »  VJnCXpf-UeK,-  ut)t\b 

+  V;sexp|i(«sc-«)ilf.  (la) 

i b-  V*  exp[i(eba  -  «)/ja‘  +  expire*  -  e)f Jr. 

(lb) 

»c»  V;hexpf-i(6*c-«)Mft 

+  V’ac  exp(-i(e*  -  e)t]a .  (lc) 

Here  the  rotating  wave  approximation  has  been 
introduced,  w  denotes  the  laser  frequency,  tie** 
is  the  excitation  energy  of  the  nucleus  (Z.  S). 
*«*■  M(Z.  N)  -  M(Z  ±  1.  JV  *  1)  is  the 
difference  of  the  nuclear  masses  and  e*  = 
f„  +  Cba.  In  the  case  of  beta'-decay  he  denotes 
the  sum  of  the  energies  of  the  outgoing  neutrino 
and  electron  or  positron,  including  the  rest 
mass.  For  orbital  electron  capture  he  includes 
the  energy  of  the  emitted  neutrino  and  the 
atomic  binding  energy  difference  of  the  cap¬ 
tured  electron  minus  the  electron  rest  mass.  At 
the  intensities  we  need  for  the  examples  at  the 
end  of  this  letter  we  assume  that  the  nucleus  is 
not  completely  ionized  so  that  electron  capture 
can  still  take  place.  We  also  assume  that  inter¬ 
nal  conversion  from  b  to  a  or  the  inverse 
process  of  nuclear  excitation  by  electron  tran¬ 
sition  [4]  are  negligible  with  respect  to  the  sti¬ 
mulated  electromagnetic  transition.  (iVab  is  then 
a  definite  electromagnetic  multipole  matrix 
clement  or  the  sum  of  several  of  them,  which 
can  be  related  to  the  half-life  of  the  level  b. 
hVcb  and  h  Vx  are  the  weak  interaction  matrix 
elements  which  govern  the  decays  from  a  and  b 
to  c.  respectively.  These  matrix  elements  include 
not  only  nuclear  wave  functions  but  also  the 
radial  part  of  the  electron  and  neutrino  state. 

We  do  not  have  to  include  phenomenological 
dampings  for  the  levels  a  and  b,  since  our  des¬ 
cription  includes  all  possible  decays  of  them.  If 
c  is  an  excited  state  we  should,  in  principle,  add 
a  damping  term.  Since  however,  the  population 
of  this  level  will  remain  small  in  any  event,  we 
can  safely  neglect  this  small  amount  of  leakage 
out  of  the  system. 

Since  due  to  the  stimulating  laser  field  the 
coupling  between  a  and  b  can  be  quite  strong 
we  shall  proceed  by  first  solving  eqs.  (la)  and 
(lb)  exactly  after  putting  c  =  0  and  then  in¬ 
tegrating  eq.  (Ic)  after  inserting  the  solutions  for 


a  and  b.  With  the  initial  condition  a( 0)  =  1. 
h(tt)  =  0  we  get  after  the  first  step 

a  =  (l/2.\)[— A_  exp(iA./)+  A,  exp(iA.r)] ,  (2a) 

b  =  ( VatJ2\)[cxp(-iA.t)~  exp(-iA.t)] .  (2b) 


A«[(*.-«)!/4  +  |Vidl]'a.  (3) 

A.  -  (to  —  eba)/2  ±  A  .  (4) 

Inserting  eqs.  (2a),  (2b)  we  can  now  integrate 
eq.  (lc)  subject  to  the  initial  condition  c(0)  =  0. 
In  the  limit  i-»*we  obtain  for  the  transition 
probability  per  unit  time 

\c\-lt  =  (7r/2.\:) 

MlV^IWWesc-e  +  A-) 

+  fifebc-  e  +  A.)] 

+  |V'J5(A:d(eac-e- A.) 

+  Aifi^-e-A-)]}.  (5) 

We  integrate  eq.  (5)  over  the  phase  space  of  the 
neutrino  and  the  electron  or  positron,  respec¬ 
tively.  for  beta--dccay 

1 10^  =  4^  /  d6‘  de-I«*“  )*]'««.«;  • 

or  just  over  the  neutrino  phase  space  for  elec¬ 
tron  capture 


J  (2nii)'  2jtVJ  d6,e: 


with  e,.  =  EJh.  e,.  =  EJU  in  order  to  obtain  the 
total  transition  rate  per  unit  time.  We  find  for 
bcta'-processes 

F  =  8^'TV  /  )*T : 


x  il  V'J.|-[A:.(eJC-  e,.  -  A.);  +  A:.(e.„  -  et  -  A_)^] 

+  |vlh|-|vv!:[(€ha-6t+A.)-: 

+  (€*-6,+  A:.).)}.  (6a) 

and  for  electron  capture 

T  =  (4w.W,)',{| V’J:[A?(eac+  e,-  A.):. 

+  A-^c^  +  e,- A.); 

+  I  VJ-|  Vbc|-[(«bc  +  ec  +  A*)J 
+  («*  +  «,+  A.):]}. 


(6b) 


Volume  I31B.  number  1.2.3 


PHYSICS  LETTERS 


10  November  1VK3 


where  xi  =  x:  for  x  >0  and  xl  =  Q  else.  The 
electron  energy  for  capture  is  given  by  et  - 
mcz/h  -  e.  where  fteb  denotes  the  atomic  binding 
energy.  The  first  term  in  eqs.  (6a),  (6b)  describes 
the  direct  decay  from  a  to  c,  whereas  the  second 
specifies  the  decay  via  the  virtual  intermediate 
state  which  can  be  significantly  enhanced  if  the 
incident  photon  is  resonant. 

Two  limits  of  eq.  (6a)  are  of  interest:  (a)  if 
the  interaction  V*  is  weak  or,  equivalently,  far 
off  resonance,  i.e..  |  V.J  <  -  »|.  we  have  for 

beta'-decay  (analogous  equations  hold  for  elec¬ 
tron  capture) 

r  =  TPP  \  (mcI/fi)J],/J 

xflvjH#,-*.)* 

+  [|VA|J|VlbP/(«b.-«)2] 

X  I(«bc  -  *')i  +  (<«  +  W  -  €«);]} .  (7) 

where  the  first  term  specifies  the  decay  rate  in 
the  absence  of  the  laser;  (b)  in  the  opposite  case 
of  comparatively  small  detuning,  i.e.  u>|«* 
V,  we  have 

r  =  HiPc*  I  d<<<*l**  ~ 

x{|V«P[(««-€.  +  |Vj)J 
+  («--«.- 1  VJ)*J 
+  |VbcP((«t,--ee  +  |VJ)i 
+  (eb,-«.-|Vlb|)?]}.  (8) 

In  this  case,  the  two  level  system  (a,  b)  is 
saturated.  Hence  the  levels  a  and  b  are  ap¬ 
proximately  equally  populated,  the  matrix  ele¬ 
ment  V*  has  cancelled  from  eq.  (S).  and  the 
direct  decay  from  a  to  c  is  reduced  by  a  factor 
of  -2  with  respect  to  eq.  (7)  (note  that  the 
term  |  V^l  in  the  square  bracket  will  be  negligi¬ 
ble  under  almost  all  circumstances!. 

Obviously,  case  (b)  is  of  paramount  interest 
since  it  allows  for  a  comparison  of  the  weak 
matrix  elements  V*  and  VN.  The  branching 
ratio  of  the  two  processes  a  -» b  —  c  and  a-»c  is 
proportional  to  |  V^p/I  V/1C|:.  The  decay  enhance¬ 
ment  due  to  the  laser  reads  in  the  particular 
case  where  V*  is  classified  as  forbidden  where¬ 
as  Vik  is  allowed, 


d  =  iKcl1  (phase  spaceW 
*“  2\  VJ:  (phase  space)*  ‘  ' 

In  addition  to  the  first  factor  in  this  ratio,  also 
the  second  can  be  much  larger  than  one.  if 
<bc  =*  +  w  >  «*.  For  allowed  beta-decays  the 

phase  space  is  given  by  the  shape  factor  (hea 
denotes  the  maximum  released  energy)  [5] 

•0 

/(£„)=  f  dep(Z.  R) 


x  «(«-’- (me  VA)2)'^*  o~ef. 

p(Z.  R)  is  the  electron  density  at  the  nucleus 
which,  if  relativistic  Coulomb  corrections  are 
taken  into  account,  depends  also  on  the  nuclear 
charge  and  radius. 

Case  (b)  applies  when 

KI  =  4|V*F/(«>.-w)J>>l.  (10) 

In  any  experiment  the  laser  will  be  tuned  to 
resonance  as  far  as  possible.  Hence,  we  should 
in  the  definition  of  K  replace  |*i..  -  <o)  by  a 
quantity  lot.  which  is  the  bandwidth  of  the  laser 
or  the  width  of  the  nuclear  level  b.  whichever  is 
larger.  The  matrix  element  |  VJ  is  related  to  the 
spontaneous  multipolc  matrix  element  |  V’JJ  by 

|V*|«N|VM.  (II) 

where  N  is  the  number  of  laser  photons,  and 
V2b  determines  the  spontaneous  half  life  of  the 
level  b  via  [6J 

|Vj,F  =  (8ir:«*,/„JV0-,c,ln2.  (12) 

with  V  the  normalization  volume.  Here  we 
assumed  again  that  internal  conversion  be  neg¬ 
ligible.  If  it  is  not.  the  lifetime  I, ;  should  be 
corrected  for  according  to  the  tabulated  internal 
conversion  coefficients  (6l_  Using  eqs.  (II)  and 
(12)  we  now  have 

K1  =  (2jt:)'ic’  In  2(.V/n(«-\W\:)-' .  (13a) 

or  in  terms  of  the  intensity  of  the  laser  field 

K'=  1.1  x  l<>-4>ic-7MV/W/(\V/cm:] 

x  (fi  :|nsl)‘'  ■  (13b) 

The  requirements  which  have  to  be  met  by  the 
laser  in  order  to  achieve  K:  2  1  are  very  severe. 
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However,  even  if  K2<*  1,  eq.  (7)  shows  that  the 
decay  is  still  enhanced  by  a  factor  of  1  +  R  with 

R  -  K2R^ ,  (14) 


a 


(68%)  (32%) 
EC  B~ 


0‘ - 071 

O* - 00 

138  a  136,  n  138  c 

56 H0  57 La 

Fiji.  2.  (a)  Level  scheme  for  the  beta-decay  'ill  —  'i!Xc. 
F.nercv  level  assignments  arc  in  McV.  Data  are  taken  from 
ref  |  ,  ,h)  Same  as  (a)  for  electron  capture  from  '!-La  to 
'sBa.  (e)  1  *-l  .ii  decays  via  a  second-order  forbidden  dccav 
both  by  /)'  and  bv  EC.  Both  decays  are  allowed  from  the 
(.1)'  level 


which  can  be  of  order  unity  or  larger  even  for 

K2<  1. 

We  now  want  to  discuss  some  examples  for 
the  scheme  discussed  thus  far.  'SI  decays  by  a 
highly  forbidden  decay  to  a  low-lying  excited 
state  of  'SXe.  From  its  excited  5/2* -state  an 
allowed  decay  to  the  same  level  can  take  place 
as  well  as  a  forbidden  decay  to  the  ground  state 
of  Xe  (see  fig.  2a).  If  we  take  for  the  parameter 
Ao>  the  natural  line  width  of  the  5/2*-level  and 
have  in  mind  that  the  half  life  f|fl  to  be  used  in 
eq.  (12)  exceeds  the  actual  half  life  by  the  in¬ 
ternal  conversion  ratio,  we  find  from  eq.  (13b) 
that  we  have  K2  =  t  for  /  =  2  x  10U  W/cm2.  A 
similar  example  relating  to  electronic  capture  is 
given  in  fig.  2b.  With  analogous  assumptions  we 
need  in  this  case  /  =  6  x  lO1’  W/cm2  in  order  to 
get  K2  =  1.  Already  a  preliminary  search  of 
nuclear  data  easily  yields  examples  with  even 
lower  ratios  hotline2.  However,  this  advantage 
tends  to  be  more  than  compensated  by  the 
rapid  increase  in  the  internal  conversion  ratio 
for  decreasing  at.  Finally,  fig.  2c  exhibits  an 
example  in  which  a  highly  forbidden  transition 
is  rendered  allowed  by  E2-absorption  of  one 
photon  with  hat  =  0.072  McV.  In  this  case,  an 
intensity  of  4  x  low  W/cm2  is  required  in  order 
to  achieve  K2  =  1.  In  all  the  preceding  examples 
the  degree  of  forbiddenness  is  reduced  by  two. 
at  least.  Hence  a  significant  enhancement  would 
already  show  up  for  K2<  1.  Whereas  the 
required  intensities  suggested  by  these  examples 
do  not  seem  to  be  extraordinary  the  assumed 
degree  of  monochromaticity  might.  Our 
assumption  of  neglecting  the  laser  band  width 
with  respect  to  the  natural  line  width  of  the 
nuclear  intermediate  level  comes  up  to 
demanding  Aa>/w  ~  8x  10“**  for  '??La  (fig.  2c) 
and  A utiut  -  1.4  x  Hr12  for  '$1  (fig.  2a).  Of 
course,  in  view  of  eq.  (13b)  the  requirements  as 
to  Aw  u>  can  be  relaxed  at  the  expense  of  a 
higher  intensity. 
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It  is  shown  that  in  the  quasiclassical  regime  total  decay  rates  of  neutral  particles  are  unaffected  by  the  application  of  an 
intense  laser  field  and  decays  of  otherwise  stable  charged  particles  cannot  be  induced.  The  method  can  also  be  applied  to 
scattering  problems  in  the  presence  of  a  laser  field. 


It  is  well  known  [  1  ]  that  in  the  presence  of  a  re*  has  been  somewhat  puzzling.  In  this  note  we  shall 

coiling  agent  there  can  be  a  significant  energy  transfer  point  out  the  reason  for  the  fact  that,  while  decay 

between  a  charged  particle  and  a  laser  field.  This  is  in*  spectra  and  differential  cross  sections  easily  exhibit 

dependent  of  the  frequency  of  the  laser  field  and  pro*  enormous  effects  due  to  an  applied  laser  field,  total 

portional  to  the  square  root  of  the  intensity.  For  an  decay  rates  and  cross  sections  stubbornly  ignore  its 

intensity  of  10 18  W/cm2  its  maximum  is  of  the  order  presence  up  to  an  intensity  of  the  order  of  the  critical 
of  1  MeV.  While  such  an  intensity  is  certainly  high,  it  one.  The  reason  will  be  found  in  the  basic  classical  ty 

is  not  out  of  reach  of  present  day  lasers  12].  Hence  of  the  interaction  with  the  laser  field  as  well  as  gauge 

the  idea  was  near  at  hand  to  investigate  whether  decays  invariance.  To  be  specific,  we  shall  concentrate  on  a 

of  elementary  particles  which  proceed  slowly  or  not  plane  wave  laser  field  of  arbitrary  polarization  and 

at  all  due  to  lack  of  energy  could  be  enhanced  or  in-  frequency  decomposition  so  that  the  explicit  Volkov 

duced  by  the  application  of  an  intense  laser  field.  solution  ( 1  ]  is  available.  However,  our  results  are 

However,  early  attempts  |3]  (for  a  recent  review  of  certainly  valid  for  a  much  wider  class  of  external 

various  decay  processes  in  various  external  fields  see  fields  including,  e.g.,  constant  magnetic  fields.  Though 

ref.  [4|)  to  calculate  these  effects  brought  about  our  approach  will  be  completely  relativistic  we  shall 

significant  effects  on  decay  rates  only  for  extremely  make  explicit  reference  to  nonrelativistic  potential 

intense  fields  of  the  order  of  the  critical  field  strength  scattering  in  the  Born  approximation  for  the  sake  of 

Eq  =  =  1.3  X  10 16  V/cm  corresponding  to  simplicity.  A  more  detailed  account  of  the  method 

an  intensity  of  4.7  X  1029  W/cm2.  The  required  in-  will  be  given  elsewhere. 

tensity  can  reduce  by  several  orders  of  magnitude  if  The  Volkov  solution  of  the  Dirac  or  Klein-Gordon 

very  small  mass  differences  are  involved  (3,4],  but  equation  in  the  presence  of  an  external  field,  speci- 

still  remains  much  larger  than  the  nowadays  attainable  fied  by  the  vector  potential  A^{x)  with  b^A^  -  0, 

1018  W/cm2.  In  view  of  the  above  mentioned  large  can  be  written  as  (we  use  four-vector  notation  with  a 

energy  transfer  at  such  an  intensity  this  non-effect  metric  such  that  ab  =  aQb0  -  ab\ we  let  c  =  1  but 

keep  h ) 
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„  t 

yp(t)  + *A  2«  )i  •  (2) 

i 

Here  p  is  the  momentum  of  the  particle  outside  the 
field,  p2  =  m2,  n  is  the  propagation  vector  of  the  laser 
field,  n2~nA=  0,  and  6  =  nx.  The  lower  limit  of  the 
integral  in  eq.  (1)  is  immaterial  because  it  merely  in* 
traduces  a  phase  which  cancels  from  all  matrix  ele¬ 
ments.  Dp(x)  is  a  Dirac  spinor  whose  explicit  form 
need  not  concern  us  here  and  which,  of  course,  is  ab¬ 
sent  in  case  of  the  Klein-Gordon  equation. 

Consider  now  a  decay  process  such  that  an  initial 
particle  described  by  the  wave  function  il/p(x)  decays 
into  N  final  particles  with  wave  functions  ^pR(x)  («  * 
1, ..  JV).  The  wave  functions  are  given  by  eq.  (1 )  or  are 
just  plane  waves  for  charged  or  neutral  particles,  re¬ 
spectively.  The  corresponding  matrix  element  is 


7  =  j  d uBfujeWl")/*. 


(7) 


“o 

In  the  quasiclassical  limit  ft  -*  0  the  stationary  phase 
method  yields 


7*  Etxp[iA(un)lti]B(un)[2mh/A"(un)}U2  ,  (8) 

n 

where  the  sum  is  over  all  zeroes  uR  of  i4*(u)  within  the 
range  u0<un<  uQ.  We  shall  actually  need  lfl2.  In 
squaring  eq.  (8),  we  can  discard  the  nondiagonal 
terms:  since  the  exponent  in  eq.  (8)  is  inversely  pro¬ 
portional  to  ft,  already  small  variations  in  the  momenta 
and  the  laser  field  strength  give  rise  to  rapid  fluctua¬ 
tions  of  the  exponential.  Hence  due  to  unavoidable 
uncertainties  in  these  parameters  the  nondiagonal 
terms  are  wiped  out  and  we  can  write 


tp ,  (*) ...  1>rN  (x)r^p(x) ,  (3) 

where  T  specifies  the  coupling.  The  integral  is  conve¬ 
niently  performed  in  terms  of  light-like  coordinates  [5] 

u*i/V5  =  (xO-x3)/V 2,  u*(x°+*3)/n/2,  (4) 

where  we  assumed  that  the  laser  field  propagates  in 
the  x3-direction.  The  corresponding  momenta  are 

Pu  3  (P°  +  P3)/\/3 ,  Pu  =  (P°  -  P3)/V 2 ,  (5) 

so  that 

2 

PX  =  UPu  +  VPv  —  S  PjXf  . 

We  then  obtain 
M=(2rrft)36<3>(p-  Ep«) 

X  j  dn  exp^-i^pu  -T)p„Mj  u/ftj 

X  exp[i(Kp({)-  Lvpn(S)UlDpr..DpjvrDp  , 

L  '  J  (6) 

where  the  three-dimensional  5-function  contains  the 
momentum  components  pv  and  pi  (t  =  1,  2).  We  shall 
evaluate  the  remaining  integral  by  the  stationary  phase 
method.  Let  us  abbreviate  the  integral  by 


t/12  =  £|B(u,)|22irftA4>„) 

n 

=  2tr ft  J du  1P(m)|26(/|'(«))  .  (9) 

In  the  last  step  we  have  made  use  of  the  representation 
8(f(x))  =  E5(x  -  xn)Hf(x„)\ , 

where  the  sum  extends  over  all  zeroes  of/(x). 

An  obvious  example  where  these  approximations 
do  not  apply  is  the  emission  of  a  photon  by  an  elec¬ 
tron  in  the  presence  of  a  laser  field,  i.e.  high-intensity 
Compton  scattering  [6j.  In  this  case  it  can  be  shown 
that  the  quantity  A(u)  becomes  proportional  to  the 
momentum  of  the  emitted  photon.  Hence  ft  cancels 
from  the  exponent  in  eq.  (7),  and  the  condition  that 
A(u)/h  varies  rapidly  as  a  function  of  its  parameters  is 
no  longer  satisfied.  This  is  related  to  the  vanishing 
mass  of  the  photon.  Consequently,  we  have  to  ex¬ 
clude  the  emission  of  a  single  zero-mass  particle  from 
our  considerations. 

Eq.  (9)  shows  that  we  can  represent  the  decay  rate 
IA/|2  as  an  integral  over  an  “instantaneous"  decay  rate 
specified  by  a  definite  value  of  the  phase  u.  This  is  a 
consequence  of  our  quasiclassical  approximation.  This 
instantaneous  decay  rate,  though  for  an  optical  laser 
field  hardly  accessible  by  experiment,  is  certainly  a 
"truly  physical  quantity”  in  the  sense  of  ref.  [7J.  As 
such,  it  must  depend  on  the  laser  field  only  via  the 
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kinetical  momenta p^u)  -p  —  eA(y)-  This  can  be 
explicitly  shown  for  eq,  (6)  by  going  over  to  the  light¬ 
like  momentum  components^)  and,  if  necessary,  to 
a  lightlike  Dirac  algebra  [5],  If  we  now  calculate  the 
total  decay  rate  by  integrating  eq.  (6)  over  the  final 
momenta  the  laser  field  is  completely  eliminated  from 
the  final  state  by  changing  the  integration  variables 
from  pn  to  Pton,,.  Obviously,  this  elimination  relies 
on  the  representation  of  the  decay  rate  as  an  integral 
over  an  instantaneous  decay  rate. 

We  now  point  out  some  immediate  consequences 
which  are,  of  course,  restricted  to  the  regime  where 
our  approximations  apply:  (i)  the  total  decay  rate  of 
a  neutral  particle  is  unaffected  by  a  laser  field;  (ii)  if 
the  decaying  particle  is  charged,  the  total  decay  rate 
does  depend  on  the  laser  field  via  its  kinetical  mo¬ 
mentum.  If  we  would  sum  over  the  initial  momentum, 
too,  the  laser  field  would  again  be  completely  elimi¬ 
nated.  This  implies:  (a)  if  the  particle  is  stable  in  vac¬ 
uum  it  remains  stable  in  a  laser  field  (an  apparent 
counterexample,  the  decay  e  -*  ey,  which  does  occur 
in  a  laser  field  but  not  in  the  vacuun,  has  been  dis¬ 
cussed  above);  (b)  if  it  is  unstable,  the  decay  might  be 
enhanced  for  particular  momenta  at  the  expense  of  be¬ 
ing  suppressed  for  others,  so  that  the  decay  rate  integrat¬ 
ed  over  all  momenta  agrees  with  its  vacuum  value.  For 
\p\  >  e\A\,  however,  this  effect  is  expected  to  be 
small;  (iii)  for  a  laser  field  we  have/ly  =  (A°-A*)l 
y/2  =  0,  hence  Pu  kjn  =p„.  Consequently,  integration 
over  n„  d-p,„  suffices  to  eliminate  the  laser  field 
from  the  final  state.  For  a  neutral  decay  this  implies, 
that  while  decay  spectra  plotted  versus  the  energies 
p%  of  the  decay  products  exhibit  dramatic  effects  due 
to  the  laser  field,  spectra  plotted  versus  pn  v  are  un¬ 
affected;  (iv)  the  decay  of  a  neutral  particle  with 
charged  constituents  is  outside  of  the  scope  of  our 
present  considerations.  E.g.,  our  results  do  not  apply 
to  multiphotoionization  of  atoms,  which  does  occur 
at  impressive  rates  and  energy  transfers  [SJ .  Also, 
whereas  the  enhancement  of  0-decay  of  a  neutron  is 
ruled  out  (conclusions  to  the  contrary  |9)  originate 
from  insufficient  numerical  estimates),  the  enhance¬ 
ment  of  nuclear  0-decay  is  still  viable,  though  not  via 
a  final  state  interaction  [10,11]. 

It  should  be  obvious  from  our  derivation  that  it 
applies  to  a  more  general  situation  than  the  plane 
wave  laser  field  which  we  considered.  We  expect  it  to 
apply  whenever  a  WKB-approximation  for  the  wave 


function, 

tp(x)~-e(~iE,  +  iS)lft  -  00) 

is  appropriate.  If  the  classical  action  S  depends  on 
more  than  one  coordinate  the  stationary  phase  ap¬ 
proach  formulated  in  eqs.  (7)  to  (9)  should  be  re¬ 
considered.  The  limits  of  applicability  of  the  latter  are 
difficult  to  assess.  It  certainly  breaks  down  when  the 
field  strength  becomes  comparable  with  the  critical 
field  strength  Eq  which  is  a  genuine  quantummechani- 
cal  quantity.  In  that  event  the  quasidassical  approach 
can  no  longer  be  expected  to  be  sufficient.  However, 
as  mentioned  at  the  beginning,  this  is  far  beyond  the 
realm  of  laboratory  experiments  which  are  within 
reach  in  the  foreseeable  future. 

There  is  a  close  connection  between  our  present 
approach  and  the  infrared  problem  of  quantum  elec¬ 
trodynamics  (which  is  classical,  too).  Cross  sections 
have  been  shown  to  be  essentially  independent  of  the 
presence  of  an  arbitrary  number  of  soft  photons  in 
the  initial  as  well  as  in  the  final  state  [12). 

Finally,  as  an  explicit  example  we  shall  consider 
non  relativistic  potential  scattering  in  Born  approxima¬ 
tion,  in  order  to  discuss  the  stationary  phase  approach 
in  more  detail  and  to  outline  a  simple  way  to  derive 
the  cross  section.  The  transition  amplitude  (3)  is  in 
tlus  case 

F= -l  fd*x  rp.{x)V{x)  ^p{x)  ,  (11) 

where  we  now  use  the  nonrelativistic  limit  of  the 
Volkov  solution  (1),  (2)  in  the  dipole  approximation. 

tpix)  =  exp|-ipx/ft 

r 

+  dt'l2pA(c')-eAH'')\J.  0=) 

The  introduction  of  lightlike  coordinates  is  now  ob¬ 
solete.  and  eq.  (9)  yields  tire  transition  probability 

\F]z  =  f  11V»-1(P-/»'))12 

T/2 

Xf  d  p0  +  (e/m)(p-p')A(t)),  03) 

-T/2 

withp0  =  p-/2m,  p'q  -  p'~l2m.  Eq.  (13)  can  be  re¬ 
written  as 
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Ifl2  =  (4>™/»)l V{h-X(p^  ~  p'^))\2 

xf  d tSip'^-pln).  ‘  (14) 

This  makes  the  exclusive  dependence  on  the  kinetical 
momenta  obvious  which  has  been  stated  above  on 
general  grounds.  We  evaluate  eq.  (13)  explicitly  for 
A 1  =  a  cos  co/,  A  2  =  -a  sin  ut.  The  result  is 

\F]2  =  (2T/h)\V(h-l(p  -  p'))\2 

x  [(eaqj/m)2  -(p'0-  P0)2)+112  ,  (15) 

where  <}T  is  the  transverse  momentum  transfer 

«t=(^(p«'-p.)2)1/2  . 

and  =  xv  if  x  >  0  and  zero  otherwise.  We  note  that 
in  the  limit  a  =  0  when  the  laser  field  is  switched  off 
we  have 

Ifl2  =  (2*77/1)1  K(/t-l(p  -  p'))|25(Po  -  P0)  ■  (16) 

The  cross  section  is 

d0 = £^4  k«*t/*>2  -  (po  -  p0)2);m 

X  |l/(/i_1(p  -  p’))l2  df2’dpo 

=  *  R  ^t/w)2  -  (p’o  -  Po)2)  ;l/2 

doel  ,  , 

X  ^  dn  dp0  ,  (17) 

which  has  to  be  compared  with  the  exact  result  [13, 
141 

In'l  n  dcjpi 

do" =  Ry"  ( eac>T/hum)  dn1  dn'  ■  0  8> 

The  latter  is  the  cross  section  under  the  condition 
that  n  photons  are  emitted,  i.e.  nf iu>  =  p0  -  p'0.  In 
view  of  eqs.  ( 1 7)  and  ( 1 8)  we  realize  that  our  ap¬ 
proximations  come  up  to 

J2(z)*n _I(r2  -  n2)~l/2  .  (19) 

Eq.  (19)  is  obtained  from  Debye’s  asymptotic  formu¬ 


las  if  (i)  n  >  1 ,  (ii)  the  actual  rapid  oscillations  of  the 

Bessel  functions  7„(c)2  as  a  function  of  n  are  averaged 

over  and  (iii)  a  rapidly  decreasing  exponential  tail  for 

n  >  z  is  neglected.  If  we  are  not  interested  in  do„  for 

a  particular  n  but  only  in  the  sum  E"  *  do„  with 

n  *  n  g i  ” 

W1  ~  n0  ^  1>  die  assumptions  (ii)  and  (iii)  are  very 
well  satisfied.  Eqs.  ( 1 7)  and  ( 1 8)  then  agree  if  we 
identify  do/dp'0  with  d o„/hu>.  The  approximations 
which  are  inherent  in  eq.  (19)  have  been  extensively 
discussed  in  ref.  [13].  Hence  we  will  be  content  with 
pointing  out  that  it  is  a  classical  approximation  ,  both 
the  rapid  oscillations  of  J„(z)2  as  a  function  of  n  and 
the  exponential  tail,  which  originates  from  energy 
transfers  rt/tto  in  excess  of  the  classical  interaction  en¬ 
ergy  eaqj/m,  are  quantum  features.  As  such  they  have 
dropped  out  of  our  quasiclassical  approximation. 

We  benefitted  from  discussions  with  M.  Hillary. 
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ABSTRACT 


It  is  shown  that  radiation  fields  well  below  the  critical  field  strength 
cannot  induce  any  noticeable  enhancements  of  nuclear  decay  rates. 
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In  a  previous  note  [1]  we  have  shown  that  the  total  decay  rate  of  a  neutral 

free  particle  decaying  into  various  charged  particles  is  unaffected  by  the 

presence  of  an  external  electromagnetic  plane  wave  field  as  long  as  the 

quasiclassical  regime  applies.  This  is,  roughly  speaking,  the  case  if  the 

external  field  is  significantly  weaker  than  the  critical  field  strength 

Ecrit  =  m2c3/eh  =  1.3x10 16  Volts/csn,  Which  is  completely  out  of  reach  for  all 

plane  wave  fields  Which  can  be  generated  in  a  laboratory.* 

Recently,  it  has  been  claimed  that  forbidden  nuclear  beta-decay  can  be  very 

significantly  enhanced  by  the  application  of  an  intense  radio  frequency  field 

with  a  field  strength  of  10  4  ...  10 5  Volts/cm  C2-4].  Hie  idea  was  that  due  to 

the  presence  of  the  external  field  the  normal  selection  rules  should  lose  their 

relevance.  Ihe  present  authors  contested  this  claim  on  the  basis  of  independent 

calculations  in  the  nonrelativistic  limit  Which  showed  that  the  enhancement 

could  be,  at  best,  proportional  to  the  ratio  E/Ecrit  of  the  applied  over  the 

critical  field  [5].  More  recent  explicit  calculations  demonstrated  that, 

actually,  the  enhancement  goes  with  (E/Ecrit)2  [6].  We  also  pointed  out  the 

analogy  with  the  above  mentioned  no-go  theorem  for  free  particles  [5].  It  has 

been  argued  that  this  theorem  since  it  was  formulated  for  free  neutral  particles 

does  not  apply  to  nuclear  decays  [7].  Tnis  criticism  is  justified,  in 

principle.  It  is  the  purpose  of  this  note  to  extend  the  no-go  theorem  to  the 

latter  case,  thus  giving  a  very  simple  and  very  general  argument  against  the 

possibility  of  altering  a  nuclear  decay  rate,  be  it  forbidden  or  not,  by  an 

external  electromagnetic  field  under  the  cited  conditions. 

The  essence  of  our  previous  argument  can  be  condensed  into  a  few  sentences: 

in  the  quasiclassical  regime  the  differential  decay  rate  can  be  written  as  an 

*  Actually,  what  matters  here,  is  the  field  strength  in  the  rest  frame  of  the  ' 
particle.  Hence,  the  quasiclassical  regime  does  not  apply,  for  example,  to 
ultrarelativistic  electrons  (y  >  10 3 )  in  an  intense  focussed  laser  field  with 
10 10  Volts/cm. 


integral  with  respect  to  time  over  an  instantaneous... decay  rate  which  depends  on 
the  instantaneous  value  of  the  vector  potential  A(t)  of  the  external  field. 

This  instantaneous  rate,  being  a  gauge  invariant  physical  quantity  [8]  must 
depend  on  the  external  field  only  via  the  kinetical  momenta  p  -  eA(t)  of  the 
charged  particles  in  the  final  state,  since  the  decaying  particle  was  assumed 
neutral.  In  obtaining  the  total  decay  rate,  the  integration  over  the  canonical 
momenta  p  can  be  replaced  by  integration  over  the  instantaneous  kinetical 
momenta  p  -  eA(t).  Thus  all  dependence  on  the  external  field  A(t)  is 
eliminated.  If  the  interaction  of  the  nuclei  with  the  field  is  neglected  as  it 
is  in  Refs.  2-4,  this  argument  carries  immediately  over  to  the  case  of  nuclear 
decays,  as  we  shall  now  demonstrate  in  detail.  We  shall  restrict  ourselves  to 
nuclear  beta-decay.  However,  the  generalization  to,  for  example,  internal 
conversion  or  radiative  transitions,  is  mainly  notational  and  should  be  obvious. 

Within  the  model  adopted  in  Refs.  [2  -  4]  the  decay  is  governed  by  the 
matrix  element 

r  iEft  -iert(u)r  *  _ 

M  ■v  jd4x  e  e  $f(r) 

-lE.t  ie.A(u)r 

*  gV  e  e  1  4>i  (r) 

*  iqx/h  ipx/h  * 


Vp(u)  *  exp(i/h)  Jdu’ (-2pA(u')  +  eA2(u)) 


(2) 


the  Volkov  part  of  the  electron  wave  function.  We  shall  ignore  spin  for  the 
time  being.  We  let  .c  *  1,  but  keep  h.  In  Eqs.  (1)  and  (2),  the  vector 
potential  A(u),  u  ■  t  -  z,  represents  the  external  field  which  propagates  in  the 
z-direction,  its  propagation  vector  being  n  =  (1,0, 0,1)  so  that  nA  =  0.  We  do 
not  have  to  specify  its  frequency  decomposition  nor  polarization.  The  neutrino 
momentum  is  denoted  by  q,  and  p  is  the  momentum  of  the  electron  outside  of  the 
field.  Ihe  nuclear  wave  functions  are  exp(-iE.  _t  +  ie  j£(u)r)(f..  ir)  for  the 

ifr  i,t 

initial  and  final  state,  respectively,  e.  -  being  effective  charges  [2].  The 

1 ,  r 

wave  functions  4.  -(r)  are  eigenfunctions  of  the  effective  nuclear  Hamiltonian 

1 1  * 

Ho  =  -  (h2/2M.  +  V(r)  with  eigenvalues  E.  ...  Contrary  to  appear  since ,  the 

i,r  1 9 r 

nuclear  wave  functions  do  not  incorporate  any  interaction  between  the  nucleus 
and  the  field;  the  factors  exp(ie^  ^Ar )  are  necessary  in  the  Coulomb  gauge  which 
we  have  tacitly  adopted  by  using  the  Volkov  solution  (2).  We  have  discussed 
this  point  extensively  in  Pef.  5.  Finally,  in  Eq.  (1)  the  weak  interction 
responsible  for  the  beta-decay  is  represented  by  gV. 

Changing  the  integration  variables  in' Eq.  (1)  to  r  and  u  and  rearranging 
terms  we  can  write 


M  'v* 


du  ofi(qT  +  PT  *  eA(u) ,  q.  -  q0  ♦  p.  -  p0  -  Hf  ♦ 


x  e 


(i/h)(Ef  -  E.  ♦  q0  ♦  p0)u 


yu) 


(3) 


where 


pfi(£T,  k2)  =  Jd3r  ^(rJgV^Or)  e 


T  T 


e 


z 


(4) 


and  we  have  used  that  e^  -  =  e  [2],  which  is  a  consequence  of  charge 

/ 

conservation.  The  subscript  T  in  Eqs.  (3)  and  (4)  designates  two-vectors, 
transverse  to  the  propagation  direction  of  the  field,  e.g.  pT  *  (p.#p.)* 

Because  the  integration  in  Eq.  (4)  is,  due  to  the  presence  of  the  nuclear  wave 
functions,  restricted  to  |r|  R0  with  R0  the  nuclear  radius,  the  fora  factor  (4) 
is  a  slowly  varying  function  of  its  arguments.  Hence  we  may  evaluate  Eq.  (3)  by 
stationary  phase  Integration,  treating  p^  as  slowly  varying.  With  obvious 
abbreviations  we  then  have 


where  the  sum  goes  over  all  zeroes  un  of  E' (u) .  The  following  is  strictly 
parallel  to  Ref.  1:  in  squaring  fl  we  can  drop  cross  terms  which  would  yield 
rapidly  fluctuating  contributions  which  average  out.  Consequently 


|M| 2  'v,  2Trh  \ 
n 


MV'2 

E"(un) 


=  2  Tth 


o(u)|25(E'(u))  . 


(6) 


so  that  the  differential  transition  rate  is  represented  as  an  integral  over  an 
instantaneous  rate.  As  argued  above,  the  latter  can  only  depend  on  A(u)  via  the 
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INTRODUCTION 

It  is  the  purpose  of  this  lecture  to  study  the  influence  of 
a  strong  external  electromagnetic  field  on  the  spectrum  of  the 
emitted  particles  in  a  nuclear  beta  decay  and  in  particular  on 
the  lifetime  of  the  radioactive  nucleus.  The  possibility  of 
manipulating  nuclear  lifetimes  [1]  in  a  laboratory  is  very 
exciting  and  could  bring  important  applications.  The  first 
experimental  proof  of  this  effect  is  more  than  30  years  old 
[2,3].  Experimentally  accessible  are  those  nuclear  decay  modes 
which  involve  an  interaction  between  the  nucleus  and  the  atomic 
electrons  such  as  internal  conversion  [3]  and  orbital  electron 
capture  [4]  which  is  closely  related  to  beta  decay.  By  changing 
the  chemical  environment  of  the  atom,  by  applying  high  pressure 
technology,  by  optical  excitation  with  strong  fields,  by  ioniza¬ 
tion  and  implantation,  etc.  the  electronic  structure  of  the  atom 
can  be  modified,  which  results  in  changes  of  the  nuclear  lifetime 
up  to  a  few  parts  in  10^  [4]. 

However,  it  always  seemed  hopeless  to  influence  the  majority 
of  the  nuclear  decay  modes  which  take  place  without  any  inter¬ 
action  between  the  decaying  nucleus  and  its  environment.  All 
imaginable  fields  in  a  laboratory  are  weak  compared  to  the  strong 
interaction  or  the  Coulomb  field  at  the  surface  of  a  nucleus. 

A Supported  in  part  by  the  Air  Force  Office  of  Scientific  Research 


Hence,  changes  of  these  nuclear  decay  rates  have  only  been  inves¬ 
tigated  under  conditions  which  are  of  interest  in  astrophysics. 
In  particular,  nuclear  beta  decay  was  studied  (i)  under  the 
condition  of  a  statistical  equilibrium  in  stellar  interiors  where 
excited  nuclear  states  are  thermally  populated  [5],  (ii)  for 
single  photon  absorption  by  the  emitted  electron  from  a  Planck 
spectrum  at  stellar  temperatures  of  the  order  of  10 3  °K  [6],  and 
< iii )  for  a  strong  uniform  constant  magnetic  field  [7]  as  it 
exists  on  the  surface  of  a  pulsar. 

,  Due  to  the  ongoing  recent  progress  in  the  development  of 
high  power  lasers  it  is  nowadays  also  possible  to  produce  extra¬ 
ordinarily  strong  fields  in  the  laboratory.  For  example,  in  the 
beam  of  a  iid-glass  laser,  which  produces  TW  pulses,  the  intensity 
can  be  (after  focusing  the  beam  down  to,  say,  ten  wavelengths)  of 
the  order  of  10i8  W/cm2,  corresponding  to  a  field  strength  of 
about  10 1 0  V/an.  These  experimental  facilities  suggest  the 
theoretical  treatment  of  nuclear  decays  in  the  presence  of  strong 
electromagnetic  plane  wave  fields. 

The  effects  of  intense  plane  electromagnetic  waves  on 
different  quantum  processes  were  already  investigated  two  decades 
ago  [8].  The  interest  soon  focused  on  the  decays  of  elementary 
particle  like  muons  and  pions  [9]  or  neutrons  [10,11]  under  the 
influence  of  a  monochromatic  external  field.  A  common  feature  of 
Refs.  7-9  and  11  is  the  result  that  the  influence  of  the  external 
field  on  the  total  decay  rate  depends  on  the  ratio  between  the 
field  strength  and  the  so  called  critical  field  strength 
Ec  3  m2c3/eh.  This  is  the  limit  for  the  applicability  of 
classical  electrodynamics,  beyond  which  quantum  effects  are 
dominant  [12].  This  result  seems  to  indicate  that  the  influence 
of  an  external  field  of  optical  frequency  on  the  lifetime  of  an 
elementary  particle  is  a  pure  quantum  effect  which  is  very  small 
as  long  as  E  is  small  compared, to  Ec.  For  electrons  the  critical 
field  strength  is  about  1.3xl0:°  V/cm  and  therefore  still  far  out 
of  the  range  of  present  laser  systems. 

On  the  other  hand,  Ref.  10  predicts  a  measurable  change  of 
the  lifetime  of  free  neutrons  in  presently  feasible  fields. 
Furthermore,  during  the  last  few  years  nuclear  decays  in  the 
presence  of  intense  plane  wave  fields  have  also  been  investi¬ 
gated.  These  calculations  predicted  appreciable  enhancements  of 
the  total  decay  rates  in  the  field  of  available  lasers,  both  for 
nuclear  gamma  decay  [13]  and  nuclear  beta  decay  [14],  and  also 
recently  for  forbidden  nuclear  beta  decay  [15]. 

In  these  lecture  notes  we  will  follow  Ref.  14  in  describing 
the  decay  process.  This  implies  the  use  of  a  modified  version  of 
the  Keldysh  approximation  [16]  which  was  introduced  for  the 
theoretical  description  of  laser-ionization  in  the  so  called 
electric  field  gauge.  This  approximation  includes  two  steps: 
(i)  the  interaction  between  the  bound  system  and  the  field"  is 


neglected,  i.e.  we  consider'  the  nucleus  to  be  unaffected  by  the 
external  field?  (ii)  in  the  final  state  the  interaction  between 
the  emitted  particle  and  the  residual  bound  system  is  neglected, 
i.e.  we  neglect  the  Coulomb  interaction  between  the  emitted  elec¬ 
tron  and  the  residual  nucleus.  The  quality  of  the  first  approxi¬ 
mation  will  be  demonstrated  in  the  next  Section.  With  the  second 
approximation  we  neglect  a  usually  small  effect  [17].  In  the 
framework  of  this  approximation  only  the  electron  emitted  in  the 
nuclear  beta  decay  couples  to  the  external  field.  This  is  why  we 
choose  nuclear  beta  decay  as  the  most  promising  decay  mode:  the 
smaller  the  mass  of  the  emitted  charged  particle,  the  stronger 
its  coupling  to  the  external  field  and  hence  the- stronger  the 
laser  inpact  on  the  decay  process .  . . . .  ’  _  .  _ _ 

There  cure  two  general  reasons  why  we  might  expect  a  change 
in  the  nuclear  lifetime  in  the  framework  of  this  model:  (i)  it 
is  well  known  for  ordinary  beta  decay  that  taking  into  account 
the  Coulomb  corrections  between  electron  and  nucleus  changes  not 
only  the  spectrum  of  the  electrons,  but  also  the  nuclear  lifetime 
[17].  Hence,  including  the  interaction  of  the  electron  with  the 
laser  field  could  result  in  the  same  effect?  -(ii)  as  we  will  show 
below  the  most  likely  energy  transfer  from,  say,  a  Nd-laser  with 
an  intensity  of  I  ^  10 18  W/cm2  to  an  electron  emitted  in  a 
typical  nuclear  beta  decay  is  of  the  order  of  MeVl  Such  a  large 
energy  increase  leads  to  a  much  larger  phase  space  of  the  emitted 
electrons  which  should  result  in  a  considerably  faster  decay  of 
the  nucleus.  . . 


In  this  lecture  we  shall  restrict  the  discussion  to  (i)  the 
nonrelativistic  theory  and  (ii)  to  electromagnetic  fields  of 
circular  polarization.  This  will  simplify  the  calculations  so 
much  that  they  can  easily  be  followed  in  detail.  Strictly 
speaking,  the  nonrelativistic  limit  holds  only  if  the  mechanical 
momentum  p-eX  of  the  electron  in  the  field  is  small  compared  to 
me.  This  implies  a  limiting  condition  for  the  kinetical  electron 
energy  E  outside  the  field,  E  =  p2/Zm  <<  me2  and  for  the  strength 
of  the  external  field  e|Al/mc  =  v<<  l.  This  approximation  seems 
to  be  unrealistic  for  many  nuclei  whose  energy  release  in  the 
decay  is  larger  than  the  electron  rest  mass  as  well  as  for  the 
fields  which  we  intend  to  consider.  Hence,  in  principle,  the 
Dirac  theory  is  required  for  the  description  of  the  electrons. 
Actually,  relativistic  effects  show  up  in  both  the  electron 
spectrum  and  the  total  decay  rate.  However,  all  the  essential 
physical  features  which  are  introduced  by  the  interaction  with 
the  laser  field  are  included  in  the  nonrelativistic  theory  in  a 
very  instructive  way.  Furthermore,  as  we  shall  demonstrate,  the 
change  of  the  lifetime  in  the  presence  of  an  external  field  is 
fairly  well  described  by  the  nonrelativistic  theory  for  all  real 
nuclei  and  arbitrary  field  strengths,  not  only  in  the  nonrelati¬ 
vistic  limit  jp-eA|  <<  me.  Hence,  we  will  present  here  only  the 
nonrelativistic  theory.  However,  all  the  numerical  results  shown 
are  gained  from  the  relativistic  beta  decay  theory  for  allowed 


transitions  with  Dirac  wave  functions  and  V-A  interaction.  This 
theory  will  be  published  elsewhere. 

We  will  first  derive  the  wave  functions  of  the  different 
particles  involved  in  the  nuclear  beta  decay.  This  is  based  on 
Ref.  18.  Whhat  follows  will  be  completely  independent  of  Ref. 
18.  Next  the  partial  transition  rates  will  be  calculated  in  the 
electric  field  gauge  and  in  the  radiation  gauge.  We  will  then 
evaluate  the  electron  distributions  and  finally  we  will  discuss 
the  nuclear  lifetime. 


WAVEFUNCTIONS 

For  the  calculation  of  .  the  wavef unctions  of.  the  different 
particles  involved  in  a  nuclear  beta  decay  we  have  to  recall"  some 
points  discussed  .in  Fief.  .18.  ..We  can  obtain  the  Hamiltonian  of  a 
particle  with  charge  e  interacting  with  an  external  ^.electro¬ 
magnetic  field  from  the  Hamiltonian  H  =  p2/2m  +  V(r)  bf  a 
particle  in  a  potential  V  with  the  help  of  the  substitution 

/  H  H3-  e  \L%-t)  . 

This  procedure  yields  the  Hamiltonian 

-H  ^  --  +  eil  J-  l/(v/(2.i) 

The  vector  potential  A8  (r,t)  and  t£ie_  scalar  potential  U8(r,t)  are 
related  to  the  electric  field  E(r,t)  and  the  magnetic  field 
2(r,t)  by 

ZpJ  =  -vu%ij  A  \y,ij  ^ 

-  Vx-P&tj. 

(2.2) 

The  index  g  denotes  the  gauge  freedom  in  the  potentials.  E  and  B 
remain^unchanged  if  we  transform  from  a  gauge  (A8,  U8)  to  a  gauge 
(A»  ,  U »  )  according  to 


where  X  is  an  arbitrary  function  of  r  and  t  . 


g 

Also,  the  state  vector  >,  which  is  a  solution  of  the 
Schrodinger  equation 

ik ^  /  '£3>  '  H9  i  )P>  ,  (2.4) 

depends  on  the  gauge.  The  wave  function  in  a  gauge  g  is 

transformed  to  another  gauge  g'  by  the  unitary  transformation 


->  ■  e  U^‘'A  ffe 


(2.5) 


The  Hamiltonian  has  to  be  distinguished  from  the  energy  operator 

e*-,£  (t-eAiv).  (2.6) 

This  (unperturbed)  energy  operator  is  a  physical  quantity  in  the 
sense  of  Ref.  19,  since  it  transforms  under  a  gauge  transforma¬ 
tion  like 

^  t  C  ft 


In  contrast,  the  Hamiltonian  is  an  unphysical  quantity  since  it 
transforms  like 


H 


e  ■He 


e  5*  %  • 


This  distinction  is  important  for  the  definition  of  an  unper¬ 
turbed  state  as  discussed  in  Ref.  18. 

Throughout  these  lecture  notes_  we  will  describe  the  laser 
field  by  a  uniform  electric  field  E(t),  neglecting  the  magnetic 
field.  This  long  wavelength  or  dipole  approximation  is  well 
justified  for  nuclear,  beta  decay  since  the  wavelength  of  visible 
light  is  about  5x10  3  cm,  whereas  the  radius  of  a  nucleus  is  of 
the  order  of  5x10  “3cm.  Hence,  when  calculating  matrix  elements 
of  the  nucleus-laser  interaction  the  external  field  can  be 
considered  as  constant  over  the  integration  area.  Within  this 
long  wavelength  approximation  two  gauges  are  most  frequently 
used;  the  electric  field  gauge  (E-gauge)  with 


and  the  radiation  gauge  (R-gauge)  with 

t 

A  it)  -  -  Jar  Far;  ^  U, \r)l )  =  0  .  (2.8) 

% 

Both  gauges  are  related  by  the  gauge  transformation 

t 

'Xn^.s  =  -  Afyijr-  -  r.Jdr  £<fj  .  (2.9) 
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Por  the  nuclear  wave  functions  to  be  used  for  the  calculation  of 
the  beta-decay  we  shall  adopt  a  simple  cne-particle  shell  model 
description:  the  nucleus  is  divided  into  a  "valence"  nucleon  and 
an  inert  oore  which  is  affected  neither  by  the  field  nor  by  the 
decay  but  generates  the  potential  V(r)  for  the  valence  nucleon 
which  decays.  The  valence  nucleon  has  an  effective  charge  e ^  and 
mass  M.  Both  do  not  have  to  be  specified;  we  have  just  to  recall 
that  because  of  charge  conservation  the  effective  charges  of  the 
initial  and  final  valence  nucleon  will  satisfy  the  relation 
eN'  i  ~  eN  f  =  ie  for  beta+-decay  (e  =  -|e|  denotes  the  electron 
charge ) . 

It  is  convenient  to  use  the  E-gauge.  The  nuclear  wave- 
function  4\-(r ,  t)  is  then  a  solution  of  the  Schrodinger  equation 

The  interaction  term  -eyEr  is  at  the  nuclear  surface  of  the  order 
of  10”  -  10*leV  if  we  apply  a  field  of  1010  V/cm.  This  is 

completely  negligible  compared  to  the  nuclear  binding  potential 
or  to  the  Coulomb  potential  at  the  nuclear  surface,  which  are  of 
the  order  of  MeV.  A  noticeable  modification  of  a  bound  nuclear 
state  would  generally  require  field  strengths  close  to  the 
critical  one.  We  shall  therefore  assume  that  the  initial  and 
final  nuclear  states  do  not  interact  with  the  external  field. 
This  is  the  first  part  of  the  Keldysh  approximation. 

As  discussed  in  Ref.  18,  the  noninteracting  state  is  defined 
as  an  eigenstate  of  the  unperturbed  energy  operator  (2.6),  not  as 
a  solutign  of  the  Schrodinger  Sq.  (2.4)  in  which  the  potentials 
a  and  Ua  are  set  equal  to  zero.  In  the  E-gauge  this  distinction 
does  not  play  a  role  since  the  Hamiltonian  Hh  and  the  energy 
coincide.  The  noninteracting  nuclear  state  is  then  a  solution  of 
Eq.^X2.10)  in  which  we  drop  the  nucleus-field  interaction  term 


(2.11) 


where  is  an  energy  eigenstate 

(  ^  +  ■  E*.  §^J  •  (2'12) 

If  we  want  to  express  the  noninteracting  nuclear  state  in  the 
R-gauge  we  have  to  be  cautious.  In  order  to  obtain  an  eigenstate 
of  the  energy  operator  (2.6)  with  the  same  constant  energy 
.eigenvalue  Ev>we  must  apply  the  unitary  transformation  (2.5)  and 
(2.9).  We  then  obtain  a  wavefunction  which  is  different  from  Eq. 
(2.11) 


(*,*)  -  e  e  <£>w<r; . 


(2.13) 


From  new  on  we  shall  drop  the  superscript  in  AR  and  refer  with 
the  notation  A  to  the  vector  potential  in  the  R-gauge  (2.8). 

For  the  calculation  of  the  electron  wavefunction  we  use  the 
second  part  of  the  Keldysh  approximation  and  neglect  the  Cbulomb 
interaction  of  the  emitted  electron  with  the  residual  nucleus. 
However,  the  interaction  between  the  electron  and  the  external 
field  should  be  taken  into  account  exactly.  Hence,  we  need  an 
exact  solution  of  the  Schrodinger  equation  for  a  free  particle 
(V  =  0)  in  the  external  field  A(t). 

We  notice  now  that  it  is  more  convenient  to  solve  this 
problem  in  the  R-gauge 

fe lr*  J  (2.14) 

instead  of  the  E-gauge 
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since  the  canonical  momentum  p  of  the  free  particle  is  a  constant 
of  motion  in  the  R-^gauge,  but  not  in  the  E-gauge.  It  is  well 
known  in  classical  mechanics  [20]  that  the  i-th  component  p^  of 
the  canonical  momentum  p  is  conserved  if  neither  the  vector 
potential  As  nor  the  scalar  potential  Us  depend  on  the  i-th  com¬ 
ponent  xi  of  the  space  coordinate  r.  In  particular,  if  A^  and  l]5 


are  spatially  uniform  p  is  conserved.  This  also  holds  true  in 
quantum  mechanics:  an  operator  is  a  constant  of  motion  if  it 
does  not  explicitly  depend  05  time  and  if  it  commutes  with  the 
Hamiltonian  HK  Hence,  if  A°  and  uj  do  not  depend  on  r  the 
eigenvalue  of  the  canonical  momentum  p  is  a  constant  of  motion. 
This  is  true  for  the  Hamiltonian  HR  in  the  R-gauge  (2.8)  with  the 
long  wavelength  approximation  for  the  field.  If  we  make  the 
additional  assumption  that  the  vector  potential  A(t)  is  switched 
on  and  off  initially  and  finally,  i.e.  A(t)  =  0  for  I t |  >  to,  we 
can  identify  the  conserved  canonical  momentum  with  the  momentum 
outside  of  the  field,  which  would  be  measured  by  a  spectrometer. 
Again,  this  holds  only  true  in  the  R-gauge.  On  the  other  hand 
the  canonical  momentum  is  not  a  physical  quantity  in  the  sense  of 
Ref.  19.  Hence  its  eigenvalues  are  different  in  different 
gauges.  Especially  in  the  E-gauge  (2.7),  where  the  vector  jpoten- 
tial  vanishes  and  the  scalar  potential  does  depend  on  r,  the 
canonical  momentum  p  coincides  with  the  mechanical  momentum 
tt(  t)  =  p  -  eA  (t)  and  its  eigenvalue  is  no  longer  conserved. 


Since  it  is  always  convenient  to  exploit  the  existence  of 
conserved  quantities,  we  like  to  solve  Eq.  (2.14)  instead  of  Bq. 
(2.15).  Only  in  the  R-gauge  can  we  make  the  following  ansatz  for 
the  electron  wavef unction:  4^  is  characterized  by  the  eigenvalue 
p  of  the  canonical  momentum  and  we  assume  that  it  factorizes  into 
an  unperturbed  plane  wave  and  a  function  f(t)  which  depends  only 
on  time,  since  also  the  field  depends  in  the  long  wave  length 
approximation  only  on  time 


i  (Ft  -  J/ % 


In  the  nonrelativist ic  theory  E  is  related  to  the  canonical 
momentum  by 


t 


t 
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(2.16) 


Only  in  the  absence  of  the  field  denotes  E  the  kinetic  energy. 
By  inserting  the  ansatz  into  the  Schrodinger  equation  (2.14)  and 
using  Eq.  (2.16)  one  obtains  the  equation  of  motion  of  f ( t ) 


a* 
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Integrating  this  equation  yields  the  exact  solution  of  Eq. 


where  V  denotes  the  normalization  volume.  The  lower  limit  of 
integration  in  Eq.  (2.17)  only  contributes  a  constant  phase  to 
the  wavef unction  and  is  therefore  insignificant.  The  analogous 
solution  of  the  Dirac  equation  is  known  as  the  Volkov  solution 
[21]. 

It  is  easy  to  check  that  the  wavef  unction  (2.17)  is  an 
eigenstate  of  the  operator  of  the  canonical  momentum  with  eigen¬ 
value  p.  Since  the  operators  of  the»  canonical  momentum  p  and  of 
the  mechanical  momentum^  * R  =  p  -  eA  commute  in  the  Rlong  wave¬ 
length  approximation,  ¥  £  is  also  an  eigenstate  of  if  and  of  the 
operator  of  the  kinetic  energy  in  the  R-gauge 


£*$.*&)  (2.13) 

The  nonrelativistic  Volkov  solution  T e  has  a  time  dependent 
energy  eigenvalue,  in  contrast  to  the  non-interacting  state  h. 
(2.13)  which  has  a  constant  energy  eigenvalue  E^. 

In  order  to  obtain  the  solution  of  Eq.  (2.15)  in  the  E-gauge 

we  have  to  carry  out  the  gauge  transformation  (2.5)  and  (2.9)  on 

4^  with  the  result 
6  \  v 
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where  we  used  Eq.  (2.16).  We  note  that  4e  depends  only  on  the 
eigenvalue  of  the  mechanical  momentum  ~(t)  =  p  -  eA(t).  This  is 
a  physical  quantity  and  has_  therefore  the  same  value  in  any 
gauge.  Here  we  expressed  *  by  the  gauge  dependent  canonical 
momentum  and  vector  potential,  both  given  in  the  R-gauge. 

Throughout  these  lecture  notes  we  shall  describe  the  laser 
field  by  a  circularly  polarized  monochromatic  plane  wave  with 
frequency  w,  propagating  in  e3 -direction 


(2.20) 


E(t)  =  ^ q  (  ^■1'ucO'd  4- 

The  corresponding  vector  potential  in  the  R-gauge  reads 

Ate)  ^  Ab  {&,  toiUi  -  ea  )  A0  -  Eq.  . 
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By  inserting  this  vector  potential  into  (2.17)  and  dropping 
the  phase  which  results  from  the-  lower  integration  limit  we 
obtain 


Here  we  introduced  the  azimuthal  and  polar  angles  $  and  8  of  the 
momentum  p 


4h  +  h  c 
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The  field  dependence  is  expressed  by  the  dimensionless  parameter 
v*  which  was  introduced  earlier 
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Er  is  again  the  critical  field  strength  Ec  =  m2c3/efi  =  1.3xl016 
V/cm.  For  practical  calculations  it  can  be  convenient  to  express 
v  either  by  the  intensity  I  and  wavelength  A  of  the  field 


=  ?,  10  A  [ yk'ryX  1  Tfb/fcy*)*  ] 


or  by  the  photon  density  o  and  the  wavelength  X 
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Finally,  the  dimensionless  amplitude  z  of  the  oscillating  term  in 
the  exponent  of  Eq.  (2.21)  is  given  by 


=  v££t- 
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For  circular  polarization  the  A2-term  in  Eq.  (2.17)  is  constant, 
A^(t)  _=  Aq  ,  and  gives  rise  to  the  field  dependent  energy  shift 


v*  mc7'2  in  7 1  (2.21).  In  a  relativistic  context,  this  term  is 


usually  interpreted  as  a  contribution  to  an  effective  mass, 
becomes  apparent  from  the  relativistic  Hamiltonian 


This 
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Hence,  in  the  following,  we  shall  refer  to  this  term  as  the 
effective  mass  correction.  It  will  play  a  crucial  role  in  the 
following  Sections. 


The  p  A  term  of  Eq.  (2.17)  is  the  origin  of  the  oscillating 
exponent  exp(iz  sin)  in  7*  (2.21)  and  responsible  for  the  fact 
that  the  kinetic  energy  (2.18)  is  not  a  constant  of  motion.  It 
can  be  interpreted  in  terms  of  the  electron  undergoing  field 
induced  multiphoton  transitions.  We  see  this  with  the  help  of 
the  generating  function  of  the  Bessel  function  Jn(z) 
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This  enables  us  to  rewrite  the  nonrelativistic  Volkov  wave- 
function  in  the  form 
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Although  we  only  deal  with  classical  fields  we  can  interpret  the 
wavefunction  (2.26)  in  terms  of  an  infinite  sum  over  n  absorbed 
or  emitted  photons. 


R 

The  time  dependent  part  of  "-0  includes  an  effective  energy 
which  changes  in  discrete  steps  of  fito.  This  represents  multi- 
photon  processes  of  any  order  n.  'We  can  estimate  the  most  likely 


energy  transfer  from  the  field  to  the  electron  by  considering  the 
behavior  of  J^(z)  as  a  function  of  the  order  n  for  fixed  z.  When 
the  rapid  oscillations  are  averaged  over,  J^(z)  increases  slowly 
with  increasing  |n|  until  it  reaches  its  maximum  around  |nl  =  z. 
For  |n|  >  z,Jn(z)  as  a  function  of  n  decreases  rapidly  to  zero. 
Hence  it  is  most  likely  that  about  z  photons  are  absorbed  by  the 
electron,  corresponding  to  an  energy  transfer  zhu  =  vcpj.  In  an 
external  field  with  v  of  the  order  of  unity  the  electron  can 
therefore  absorb  energy  from  the  field  up  to  the  order  of  its  own 
kinetic  energy,  typically  seme  MeV  for  nuclear  beta  decay.  This 
strong  energy  absorption  encourages  us  to  expect  a  considerable 
laser  effect  on  nuclear  beta  decay.  However,  we  realize  already 
that  this  will  be  counteracted  by  the  increase  of  the  electron's 
effective  mass. 

Finally,  we  need  the  neutrino  wavefunction.  Since  the 
neutrino  neither  interacts  with  the  laser  field  nor  with  the 
residual  nucleus,  it  is  well  described  by  a  plane  wave  solution 
of  the  Klein-Gordon  equation 

*  17  *  <*</>[-%  (?vt- ft)  1  (22?) 

with  the  relativistic  dispersion 


TRANSITION  RATS 

We  now  intend  to  calculate  in  the  framework  of  the  model 
derived  in  the  last  Section  the  effects  of  a  strong  laser  field 
on  measurable  quantities  like  the  spectrum  of  the  emitted  elec¬ 
trons  and  the  lifetime  of  the  nucleus.  The  transition  amplitude 
reads  in  first  order  of  perturbation  theory  with  respect  to  the 
weak  interaction 


T/t 


The  weak  interaction  is  denoted  by  (gV)  where  V  is  a  dimension¬ 
less  operator  and  g,  is  the  coupling  constant  for  beta  decay, 
g  =  1.4x10  3  erg  cm".  According  to  the  Keldysh  approximation 

the  electron  wavefunction  '+'e  is  given  by  the  exact  solution 
(2.21)  or  (2.18)  for  a  free  particle  in  a  circularly  polarized 
field  and  the  nuclear  wavefunctions  L  and  Aj  p  are  given  by 
the  noninteracting  solution  (2.11)  or  (£.13),  respectively,  ihe 


uncharged  neutrino  is  described  by  the  plane  wave  (2.27).  By 
taking  the  wavefunctions  of  all  charged  particles  consistently  in 
the  same  gauge,  either  the  E  or  the  R-gauge,  we  obtain  in  any 
case  o 
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If  we  choose  the  wavefunctions  in  the  R-gauge  we  have  to  make  use 
of  the  relation  evf”e\jis_e  in  order  to  obtain  Eq.  (3.2). 
For  the  energy  balance  orie  has  to  take  into  consideration  the 
mass  of  the  electron  which  is  created  during  the  beta  decay. 
Hence  the  value  Q  in  Eq.  (3.2)  is  given  by 
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Since  we  are  only  interested^in  allowed  nuclear  beta  decay  we  can 
drop  the  factor  exp{  i  (p-eA-hq)  r/T> }  in  the  matrix  element.  The 
physical  justification  for  this  approximation  is  the  fact  the 
electron  and  the  neutrino  do  not  carry  away  any  angular  momentum 
in  an  allowed  beta  decay  so  that  the  spatial  dependence  of  their 
wavefunctions  can  be  neglected.  In  the  nonrelativistic  theory 
this  is  a  fair  approximation.  The  nonrelativistic  limit  implies 
Ip  -  eA(t) I  <<  me,  as  discussed  above.  Since  the  spatial 
integrtion  in  the  matrix  element  extends  over  the  radius  R  of 
the  nucleus,  which  is  of  the  order  of  some  10  ljcm,  the 
exponential  factor  can  be  estimated  for  all  nuclei  by 


|(j£~e4  )r/%  I  -e 4  /  fc/t  <£ 


Since  we  now  dropped  the  factor  exp  ieAr/h,  the  discussion  of  the 
preceding  Section  on  the  question  of  when  to  take  which  wave 
function,  viz.  (2.11)  or  (2.13),  appears  obsolete.  Nevertheless, 
we  have  emphasized  this  point  because  it  is  important  in 
principle  and  of  crucial  significance  in  the  case  of  forbidden 
beta-decay  [16,  32]. 

By  the  same  arguments  the  space  dependent  factor  exp(iqr/h) 
in  the  neutrino  wavefunction  can  be  replaced  by  unity.  If  the 
nuclear  matrix  element  is  abbreviated  by 


the  transition  amplitude  for  allowed  beta  decay  takes  the  form 
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With  the  help  of  the  expansion  (2.25)  Eq.  (3.4)  can  be  rewritten 
in  the  form 
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This  representation  again  allows  for  a  simple  interpretation  in 
terms  of  n  photons  emitted  or  absorbed  by  the  electron. 

For  the  calculation  of  the  decay  rate  we  only  need  the 
transition  amplitude  (3.5)  in  the  limit  T ->•  <*>: 


4-w,  F(TJ  -  -,77  1.  •JibJe 

T->  &o  h  =  -  ^ 


—  xoo  i 


trp  -  Q  ~ 


(3.6) 


With  the  help  of  the  standard  relation  [22] 
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we  obtain  the  transition  rate  per  unit  time 
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Only  diagonal  terms  n  =  n'  contribute  to  I F(T)  1 2  in  the  limit 
T  -*•  <*>  since  the  off-diagonal  terms  yield  products  of  two  delta 
functions  with  different  arguments.  Due  to  reasons  which  will 
become  clear  below  it  is  convenient  to  introduce  the  dimension¬ 
less  quantity 

■^o  =  ( *  +  3  -  & )  Abu) ) . 


Eq.  (3.7)  can  then  be  rewritten  in  the  form 
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In  the  next  step  we  will  simplify  the  sum  over  the  Bessel 
functions.  In  high  power  laser  fields  the  argument  z  of  the 
Bessel  function  can  achieve  values  of  10 3  to  10°  as  discussed 
above.  Hence,  up  to  10 5  terms  can  contribute  to  the  sum  (3.8). 
An  excellent  approximation  for  Bessel  functions  of  large  order  n 
and  for  z  >  |n|  is  Debye's  asymptotic  expansion  [23]: 

~  ji  (2  ~n  )  ^ Cos  (Jt- k  ‘  -  yi ctrccos -  \ 

The  rapidly  decreasing  exponential  tail  of  Jn(z)  for  |n|  >  z  will 
be  neglected,  i.e.  we  set  Jn(z)  =  0  for  |n|  >  z.  This  approxima¬ 
tion  is  only  exact  in  the  limit  z .  According  to  Eq.  (2.24) 
this  corresponds  to  the  classical  limit  h o.  As  long  as  z  is  a 
large  number  this  classical  limit  is  well  justified.  Further¬ 
more,  after  squaring  Jn(z)  we  can  average  over  the  oscillations 
in  (z)  since  it  is  hardly  possible  in  any  experiment  to 
determine  the  intensity  parameter  v  and  therefore  also  z  to  such 
an  accuracy  that  the  phase  of  the  Bessel  functions  is  precisely 
defined.  Also  these  rapid  oscillations  of  J  £  (z)  as  a  function 
of  n  are  a  quantum  feature  [24].  This  procedure  then  yields  a 
powerful  approximation  for  the  square  of  Bessel  functions 


The  0-function  denotes  the  usual  step  function  0(x)  =  1  for  x  >  0 
and  =  0  for  x  <  0. 

Since  we  are  not.  interested  in  a  particular  term  of  Eq. 
(3.8)  with  a  fixed  number  n  of  transferred  photons,  but  only  in 
the  total  effect  of  very  many  multiphoton  terms,  we  car*  replace 
the  sum  in  Eq.  (3.8)  by  an  integral.  By  inserting  the  approxima¬ 
tion  (3.9)  into  Eq.  (3.8)  and  carrying  out  the  n-integration  we 
obtain  the  transition  rate  per  unit  time  in  the  quasiclass ical 
limit 

If  we  recall  that  z  and  n  cure  proportional  to  (hw)  1  we  notice 
that  in  the  quasiclassical  approximation,  which  applies  in  the 
limit  hjj  <<  mc:,  the  transition  rate  w  becomes  independent  of  the 
field  frequency  u.  It  depends  on  the  external  field  only  via  the 
intensity  parameter  v  (2.23). 


ALTERNATIVE  DERIVATIONS  OF  THE  QUASICLASSICAL  TRANSITION  RATE 

There  are  two  further  instructive  and  easy  methods  to  derive 
Eq.  (3.10)  which  illustrate  its  classical  character.  In  the 
first  part  of  the  current  Section  we  calculate  the  time  integral 
for  the  transition  amplitude  (3.4)  with  the  help  of  the 
stationary  phase  method  [25]  instead  of  expanding  the  integrand 
in  terms  of  Bessel  functions.  The  method  of  stationary  phase 
approximates  the  integral  over  a  rapidly  oscillating  function  by 
the  contributions  from  the  regions  around  stationary  points  to 
the  integrand  [26]: 


Here  x  is  a  large  positive  variable  and  f(t)  a  real  function  of 
the  real  variable  t,  so  that  the  integrand  is  rapidly  oscillating 
unless  f(t)  is  stationary.  Hence  the  major  contribution  to  the 


value  of  the  integral  arises  from  the  vicinity  of  the  points  tn 
at  which  f'(tn)  =  0.  The  sum  over  n  extends  over  all  stationary 
points  in  the  interval  [a,b],  The  stationary  phase  method 
becomes  exact  in  the  limit  x  -*08 ,  i.e.  for  infinitely  large 
exponents.  Hence  the  use  of  this  approximation  method  for  the 
time  integration  in.  Eq.  (3.4)  implies  the  transition  to  the 
classical  limit  h  V  o  (i.e.  z  -*■  00 )  already  in  the  transition 
amplitude  F(t). 

In  order  to  apply  the  stationary  phase  method  to  Eq.  (3.4) 
we  first  have  to  determine  the  stationary  points  tn  of  the 
exponent  of  the  integrand: 

Cos  (urbn,  +  y  )  =  4-  (it0  +  E?  )  , 


In  order  to  obtain  stationary  points  we  find  again  the  condition 
z2  >  (n0  +  Ejhu)2.  We  furthermore  note  that  the  second  deriva¬ 
tive  of  the  exponent  has  the  same  form 
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at  all  stationary  points.  Thus  we  obtain  for  the  transition 
amplitude  in  the  semiclassical  limit 
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Cross  terms  from  different  stationary  points  will  once  again  not 
contribute  to  the  transition  rate  |F(T)|"  since  the  phases  at  the 
stationary  points  tn  are  randomly  distributed  for  all  practical 
purposes  so  that  the  cross  terms  cancel  each  other  in  the  limit 
T  -coo.  we  have  two  stationary  points  per  time  interval  2^/u  and 
hence  wT/~  stationary  points  in  the  entire  integration  region 
[-T/2,  T/2].  We  thus  find  from  counting  all  stationary  points  in 
I F(T) I  2  the  transition  rate  per  unit  time 


(4.2) 
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This  result  is  identical  with  the  quasiclassical  approximation 
(3.10)  of  the  Bessel  function  approach. 

There  is  an  even  simpler  argument  leading  to  Eq. (4.2)  which 
is  classical  from  the  outset.  We  can  write  the  transition  rate 
per  unit  time  for  the  decay  in  the  absence  of  the  laser  field  as 
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where 
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is  the  electron  energy  distribution  function  and  E  =  p2/2m  the 
kinetic  energy.  Of  course,  in  the  absence  of  the  field  the  elec¬ 
tron  energy  E  is,  for  specified  Q,  a  function  of  the  neutrino 
energy  Ev  only.  In  the  presence  of  the  field  the  kinetic  energy 
of  an  electron  is  no  longer  conserved  but  time-dependent  (cp.  Sq. 
(2.61)),  so  that 
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with  p  the  conserved  canonical  momentum.  We  now  assume  that  when 
the  field  is  switched  on,  Vfq  in  Eq.  (4.3)  is  unaffected  but  the 
electron  energy  distribution  ?free(E)  is  changed  into 
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This  is  the  time  average  over  the  instantaneous  electron  energy 
distribution.  If  we  recall  that 
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where  the  sum  is  over  all  zeroes  Xj  of  f(x),  and  notice  that 
Q-Ey-E(t)  has  two  zeroes  within  to<t  ^  to  +  2tt/w,  we  find  that 
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agrees  with  Eq.  (4.2).  This  procedure  illustrates  once  again 
that  the  approximation  (3.9)  for  the  Bessel  functions  as  well  as 
the  stationary  phase  method  correspond  to  a  purely  classical 
treatment  of  the  electron-laser  interaction. 


By  comparing  Eqs.  (4.4)  and  (4.6)  we  can  already  anticipate 
a  general  feature  of  the  electron  energy  distribution  which  will 
be  more  extensively  discussed  in  the  next  Section.  Due  to  the 
field,  the  kinetic  energy 
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in^  Eq.  (_4.4)  is  ^replaced  by  E(t)  of  Eq.  (4.5).  Now,  if 
|eA|  >>  IPfree!»  Ipi  will  be  much  larger  than  lpfreel  in  order 
that 
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can  be  satisfied  as  prescribed  by  Eqs.  (4.4)  and  (4.6).  Hence, 
for  strong  fields,  we  expect  the  electron  energy  spectrum  to 
extend  to  much  higher  energies. 


ELECTRON  DISTRIBUTIONS 

In  order  to  obtain  the  angular  distribution  and  the  energy 
spectrum  of  the  emitted  electrons  the  transition  rate  w  has  to  be 
integrated  over  the  phase  space  of  the  neutrino  and  the  electron. 
The  total  decay  rate  F  is  given  by 


It  is  convenient  to  use  polar  coordinates  for  the  momenta  p  and  q 
and  to  apply  the  dispersion  relations  (2.16)  and  (2.23)  for  the 
electron  and  the  neutrino.  This  yields 
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For  the  transition  rate  w  in  Eq.  (5.1)  we  use  the  semiclassical 
limit  (3.10)  which  was  shown  to  be  an  excellent  approximation. 
Since  w  (3.10)  depends  only  on  Ev  and  via  z  and  n0  on  E  and  6  the 
integrals  over  9  ,  <j>  and  4>  are  trivial  in  Eq.  (5.1).  After  sub¬ 
stituting  Ev  by  v  v 


We  have  two  restrictions  for  the  integration  variable  x  which 
have  to  be  satisfied  simultaneously:  (i)  n0  <  x,  corresponding  to 
the  condition  of  positive  neutrino  energy  Ev  >  0,  and  (ii) 
-z  ^  x  <  z  originating  from  the  9  -function.  Therefore  we  have 
to  distinguish  between  three  cases:  (i)  if  no  >  z  the  total 
decay  rate  r  vanishes;  (ii)  if  -z  ^  ng  <  z  the  x- integration 
extends  from  no  to  z;  (iii)  if  nn  <  -z  the  x-integration  extends 
from  -z  to  +z.  This  leads  to  the  elementary  integrals 


The  angle  and  energy  distribution  of  the  electrons  then  reads 
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Since  no  and  z  depend  on  E  and  0,  the  condition  z  >  no,  i.e. 

t r  ./-  —  -  6^  —  (!V  vM*  F  vT^'hiF 

yields  the  limits  of  allowed  emission  angles  9  and  of  possible 
electron  energies  E.  For  E  <  Q  -  \)2rac-L/2  this  condition  is 
always  satisfied.  But  for  higher  energies  electrons  are  only 
emitted  into  a  certain  angular  range.  The  limits  for  the 
emission  angle  are  for  fixed  energy  E 
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The  limits  for  the  electron  energy  can  be  obtained  from  the 
condition  sine  <  Is 
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Figure  1  shows  for  the  relativistic  theory  the  possible  (E,9) 
values  of  the  emitted  electrons  in  a  laser  field  with  an  inten¬ 
sity  parameter  v  =  0.5  for  two  examples.  18F  decays  via  B+  decay 
to  l80  with  a  moderate  energy  release  Q  =  634  keV,  3H  decays  via 
B  decay  to  3He  with  the  very  small  energy  release  Q  =  13.6  keV. 
(Note  that  we  defined  Q  as  the  mass  difference  between  nuclei, 
not  between  neutral  atoms).  The  solid  line  in  Fig.  1  shows  the 
boundary  of  the  integration  area:  on  the  left  or  inside  the 
solid  line,  respectively,  we  have  the  allowed  (E,9)  values.  On 
the  right  or  outside  this  curve  no  =  z,  respectively,  we  have  the 
nonclassical  regime  with  extremely  few  events  which  are  neglected 
in  the  present  classical  theory.  Whereas  Eq.  (5.3)  yields  a 
symmetric  electron  distribution  around  6  =90°  with  respect  to 
the  laser  beam  axis,  the  relativistic  calculation  yields  an 
asymmetric  distribution.  Electrons  at  the  high  energy  end  of  the 
spectrum  are  mainly  emitted  in  forward  direction.  The  dotted 
lines  indicate  a  value  of  no/z  >  1  up  to  which  the  Bessel 
functions  are  already  so  extremely  decreased  that  this  area  can¬ 
not  play  a  role  in  the  nonclassical  regime  anymore.  Nonclassical 
corrections  to  Eq.  (5.2)  can  only  be  due  to  (E,0)  points 
extremely  close  to  the  solid  line.  Furthermore,  the  dashed  lines 
in  Fig.  1  show  contours  of  constant  ratios  no/z<  1.  The  smaller 
this  ratio  is  at  some  point  (E,9),  the  higher  is  the  number  of 
events  that  contribute  to  d2r/dEde  at  this  point. 

For  the  example  3H  a  laser  field  with  v  =  0.5  is  already  so 
strong  that  the  effective  mass  exceeds  the  Q  value.  For  all 


Fig.  1.  Area  of  ail  lowed  emission  angle  and  energy  (solid  curve) 
and  contours  of  oonstant  ratio  no/z  for  a  field  inten¬ 
sity  v  =  0.5  and  for  the  nuclei  (a)  1SF  and  b)  3  H.  The 
relativistic  electron  energy  is  scaled  in  units  of  me3: 
e  =  (m^4  +  p2c2 )  Vmc 2.  These  curves  come  from  the 
complete  relativistic  theory. 


energies  n0  is  then  positive.  This  is  in  contrast  to  the 
field-free  situation  where  no  is  always  negative.  The  given 
combination  of  Q  and  v  allows  oily  configurations  with  n0/z>O.83. 
The  value  no  denotes  the  minimum  number  of  photons  which  the 
electron  has  to  absorb  from  the  field.  In  this  example  the  elec¬ 
tron  must  absorb  energy  from  the  field  since  the  energy  release 
of  the  nuclear  decay  alone  is  not  sufficient  to  account  for  the 
effective  mass  of  the  electron  in  the  field.  Since  n9  for  all 
(E,r)  values  is  rather  close  to  z,  only  relatively  few  multi¬ 
photon  terms  contribute  to  dzr/dEd9  at  one  particular  (E,~) 
point. 

All  these  effects  can  also  be  seen  in  the  spectrum  of  the 
electrons.  We  obtain  the  energy  distribution  of  the  electrons  by 
inte,  ?ting  Eq.  (5.2)  numerically  over  9  between  the  limits 
(5.3).  Figure  2  shows  the  results  of  the  corresponding  relativ¬ 
istic  calculation,  again  for  the  two  examples  1 3F  and  3H.  A 
similar  plot  can  be  found  in  Ref.  27.  With  increasing  field 
intensity  the  electron  distribution  extends  to  much  higher 
energies.  On  the  other  hand,  the  maximum  of  the  spectrum 


Fig.  2.  Energy  distribution  of  the  electrons  for  different  field 
intensities  for  the  nuclei  a)  loF  and  b)  3  H.  The  rela- 
tivistic  electron  energy  is  scaled  in  units  of  me"  as  m 
Fig.  1,  eq  is  related  to  the  Q  value  of  the  decay  by 
Go  =  1  +  Q/mc2,  and  Y  denotes  a  constant  factor  of 
dimension  (energy  x  time).  These  curves  are  calculated 
from  the  complete  relativistic  theory. 


decreases  with  increasing  v.  The  area  under  the  curve,  which 
represents  the  total  number  of  emitted  electrons  and  therefore 
the  lifetime  of  the  nucleus,  seems  to  be  more  or  less  constant 
for  the  different  field  intensities.  We  shall  investigate  this 
question  in  the  next  Section. 

In  the  3H  spectrumin  Fig.  2b  we  see  again  the  effect  of  the 
strong  effective  mass  v"mc"/2  exceeding  the  Q-value  of  the  reac¬ 
tion.  For  sufficiently  high  v-values  the  spectrum  does  not  begin 
at  E  =  0  anymore,  but  at  some  finite  energy.  This  reflects  again 
the  fact  that  the  electrons  from  the  low  energy  end  of  the 
laser-free  spectrum  must  absorb  a  considerable  amount  of  photons 
from  the  field  in  order  to  raise  their  effective  mass,  since  the 
available  Q  value  is  not  sufficient. 


NUCLEAR  LIFETIME 


One  possibility  to  obtain  the  total  decay  rate  r  and  the 
lifetime  of  the  nucleus  r  =  1/ F  is  the  numerical  integration  of 
Eq.  (5.2)  between  the  limits  (5.3)  for  the  angle  and  (5.4*  for 
the  energy.  However,  it  is  also  possible  to  calculate  F 
analytically  from  Eq.  (3.10)  by  choosing  another  set  of  variables 
for  the  electron  momentum  and  another  order  of  integration.  The 
dependence  of  z  on  (Pi  +  pt)  *  suggests,  for  example,  the 
following  set  of  variables 
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For  the  neutrino  momentum  q  we  use  again  polar  coordinates.  With 
this  choice  for  d  and  d3q  and  with  the  semiclassical  transition 
rate  (3.10)  the  total  decay  rate  r  (5.1)  reads  after  carrying  out 
the  three  trivial  angular  integrations 


After  expressing  z  and  no  in  terms  of  s  and  ps  the  6-function 
yields  the  following  restrictions  for  the  integration  limits  if 
the  integrations  are  carried  out  in  the  indicated  order  (first 
ds,  then  dp  3  and  finally  dEy): 
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The  total  decay  rate  "  can  then  be  rewritten  in  the  form 
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The  value  of  the  integral 
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is  independent  of  the  intensity  parameter  v  ,  and  so  are  the 


limits  of  the  remaining  integrations  over  p3  and  E, 


Hence, 


after  integration  over  ds  the  dependence  on  the  external  field 
drops  out  of  the  expression  (6.2)  for  the  total  decay  rate.  This 
means  that  in  the  quasiclass ical  limit  the  total  decay  rate  r 
is  not  affected  by  the  external  field,  nor  are  the  partial  decay 
rates  dr/dp3,  dr/d$  and  d°I/dJq.  Field  effects  only  show  up  in 
dr/ds  and  in  the  related  partial  rates  dr/dE»  dr/de,  dF/dpi, 
dr/dp2,  etc. 


The  last  two  integrations  in  Eq.  (6.2)  over  dp3  and  dEv  are 
easily  carried  out.  We  then  obtain  for  r  the  nonrelativistic 
limit  of  the  total  decay  rate  of  allowed  nuclear  beta  decay  in 
the  absence  of  an  external  field: 

- r-7  & 

o In  .—,,  J 


/ 


_  <P(/T 


Jcfeje 


(6.3) 


/CfrV/j 


(WX*  &VJ-  , 


Again,  there  is  a  simpler  way  to  realize  the  field  independence 
of  r  in  the  classical  limit:  according  to  Eqs.  (4.5-7) 


■J 


I)*/'* 

-+0+ 'Jr/co 

JcH  f  ot*£>  (P (Q  -  by  J 


(6.4) 

o? 


v 


where  the  integrations  over  t  and  p  have  been  commuted.  If  we 
now  transform  the  integration  variable  p  to  p  -  eA(t),  the  field 
dependence  is  entirely  eliminated.  Obviously,  in  order  to 
achieve  this,  it  suffices  already  to  integrate  over  the  two  com¬ 
ponents  of  p  which  lie  in  the  plane  of  A(t). 

The  relativistic  theory  also  leads  to  the  result  that  the 
nuclear  lifetime  is  unaffected  by  the  external  field.  The 
physical  importance  of  the  effective  mass  term  \j2  mc2/2  for  this 
result  should  be  stressed  again.  We  recall  that  the  effective 
mass  is  due  to  the  K2- term  of  the  interaction  Hamiltonian  in  the 
R-gauge.  If  the  field  is  strong  enough  so  that  the  effective 
ness  exceeds  the  Q  value  of  the  decay  emitted  electrons  must 
absorb  the  energy  difference  from  the  field.  We  thus  have  two 
competing  processes:  since  the  electrons  absorb  energy  they  have 
a  larger  phase  space  and  the  nucleus  should  decay  faster.  On  the 
other  hand,  due  to  the  effective  mass  oily  those  electrons  are 
emitted  which  absorb  between  ng  and  approximately  z  photons.  The 
distribution  of  the  photon  absorption  with  respect  to  the  photon 
number  n  is  given  in  the  classical  limit  by  the  approximation 
(3.9)  for  the  Bessel  functions  or  by  the  classical  distribution 
(4.6).  With  increasing  v,  i.e.  with  increasing  phase  space,  the 
minimum  of  no/z  in  the  entire  (E, 6)  plane  approaches  unity  so 
that  the  portion  (z  -  no)z  of  photon  absorptions  that  still 
leads  to  decay  becomes  very  small.  In  the  limit  of  classical 
electron  field  interaction  these  two  effects,  increasing  the 
phase  space  and  overcoming  the  effective  mass,  cancel  out  each 
other  exactly. 

In  a  different  context,  namely  the  hyperfine  splitting  of 
atomic  levels  in  a  strong  radio  frequency  field,  the  just  men¬ 
tioned  cancellation  has  been  noticed  previously:  by  merely 

replacing  the  electron  mass  m  by  the  effective  mass  (m2+e2A2/c2 ) 
in  the  Hamiltonian,  an  effect  was  predicted  which  was  not 
corroborated  by  experiment  [33].  It  was  later  shown  [34],  that 
in  the  long  wave  length  limit  this  effect  is  cancelled  by  the 
pA-term  which  was  initially  neglected.  This  is  another  example 
for  the  interplay  between  the  pA-  and  the  A2-  term.  The  argument 
following  Eq.  (6.3)  shows  clearly  that  this  cancellation  is  just 
a  consequence  of  the  minimal  coupling  interaction. 

However,  one  should  keep  in  mind  that  the  field-independence 
of  the  nuclear  lifetime  holds  only  in  the  classical  limit  ft  -*■  0 
for  the  electron- field  interaction.  This  limit  is  well  justified 
for  all  beta  decaying  nuclei  in  the  presence  of  laser  fields  with 
intensity  v  in  the  order  of  unity  and  with  optical  or  longer 
wavelengths.  However,  the  quasiclassical  limit  breaks  down  for 
fields  with  shorter  wavelengths  (x-ray  regime)  and  for  (ficti¬ 
tious)  nuclei  with  tiny  Q  values  Q  <<  mc^.  In  these  cases  the 
parameter  z  which  corresponds  to  the  maximum  number  of  absorbed 
photons  cannot  be  considered  a  large  number  anymore.  The 
classical  approximation,  however,  holds  only  for  multiphoton 


transfers  of  very  high  order.  The  classical  limit  furthermore 
breaks  down  for  super  intense  fields  when  v  >>  1.  Then  the 
effective  mass  term  becomes  so  large  that  the  minimum  value  of 
the  parameter  n0,  which  corresponds  to  the  minimum  number  of 
absorbed  photons,  comes  very  close  to  z,  i.e.  there  are  no  (E,9) 
configurations  with  no  <<  z  anymore.  In  the  Bessel  function 
procedure  this  implies  that  (E,  9)  configurations  beyond  the 
no  =  z  border  become  important.  For  the  stationary  phase  method 
the  case  v  >>  1  means  that  higher  order  terms  in  the  expansion 
around  the  stationary  points  have  to  be  taken  into  account. 
These  three  ways  to  get  away  from  the  validity  of  the  semiclassi- 
cal  limit  will  be  demonstrated  again  in  the  next  Section. 


THE  NONCLASSICAL  REGIME 


In  order  to  calculate  the  quantum  effects-- of  the  elec¬ 
tron-field  interaction  we  have  to  attempt  to  calculate  the  total 
decay  rate  analytically  as  far  as  possible  without  any  approxima¬ 
tions.  We  therefore  start  from-  the  exact  expression  (3.7)  for 
the  transition  rate  per  time  in  terms  of  Bessel  functions  and 
choose  again  the  variables  (s,$,p?)  for  the  electron  momentum  as 
in  the  last  section.  The  sum  over  the  Bessel  functions  in  Eq. 
(3.7)  cannot  be  carried  out  anymore.  We  choose  the  order  of 
integrations  so  that  we  avoid  an  integration  over  the  Bessel 
function.  The  three  integrations  over  the  angular  variables  are 
trivial,  as  well  as  the  integrations  over  the  neutrino  energy  Ev 
and  over  p  3.  We  thus  obtain  an  exact  expression  for  the  nuclear 
lifetime  in  a  laser  field 
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The  summation  over  n  and  the  integration  over  s  can  be  carried 
out  analytically.  For  this  purpose  the  sum  of  the  squared  Bessel 
functions  has  to  be  transformed  in  a  proper  way,  and  the  calcula¬ 
tion  of  the  total  decay  rate  r  is  somewhat  lengthy  [28]. 

However,  there  also  exists  a  completely  different  approach 
for  the  exact  calculation  of  the  total  decay  rate  T  in  the 
presence  of  an  external  field.  It  is  possible  to  express  the 
field  influence  on  r  by  an  operator  acting  on  the  decay  rate 
r  free  (6.3)  of  the  beta  decay  in  the  absence  of  the  field  [29]. 
.The4  exact  expression  for  the  nonrelativist ic  rate  f  of  allowed 
beta  decay  reads 


For  circularly  polarized  fields  the  operator  M  has  the  form 
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As  usual,  these  operators  are  defined  by  their  power  series 
expansion 
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The  coefficients  c^  are  given  ty 
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Noticing  that  x  occurs  in  Eq.  (7.3)  in  combination  with  u  as  the 
product  ux,  we  can  simplify  the  operator  considerably.  The  k-th 
term  of  the  operator  M  applied  to  the  free  decay  rate  rfree 
differs  by  the  order  (hu/Q) 2  from  the  (k-l)-th  term.  In  the 
limit  hcj  <<  Q  it  is  therefore  sufficient  to  take  only  the  first 
term  k  =  1  in  M  into  account.  Consistently  also  exp{-iM/h}  can 
then  be  cut  off  at  £  =  1.  To  first  order  of  (hu/Q)  the  total 
decay  rate  is  thus  given  by 
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The  enhancement  of  the  total  decay  rate  due  to  the  external  field 
is  then  to  first  order  of  (hu/Q) 


(7.5) 
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Expressing  the  intensity  parameter  v  by  the  critical  field 
strength  Ec  (2.23)  the  enhancement  can  be  written  in  the  form 
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This  result  agrees  with  the  nonrelativistic  limit  in  Ref.  28.  As 
mentioned  at  the  end  of  the  preceding  Section  we  can  obtain  a 
considerable  enhancement  either  for  nuclei  with  very  small  Q 
values  Q  <<  me2,  or  for  x-ray  fields  or  for  high  field  intensi¬ 
ties  v  >  l. 


Figure  3  shows  the  enhancement  R  (7.5)  as  a  function  of  v 
for  a  Nd  laser  with  hu  =  1.17  eV  and  for  several  Q-values.  For 
Q  <<  me2  the  enhancement  already  becomes  very  large  for  rela¬ 
tively  moderate  field  strengths.  It  should,  however,  be  men¬ 
tioned  that  the  chosen  Q  values  are  considerably  smaller  than  the 
Q  values  of  any  existing  beta  decaying  nuclei.  For  3H  Eq.  (7.5) 
requires  already  v  =  4x10 3  for  a  Nd  laser  in  order  to  obtain 
R  =  2,  corresponding  to  a  field  strength  Eo  of  1.2x10  14  V/cml 


V 


Enhancement  R  of  the  total  decay  rate  r  due  to  an 
external  field  of  wavelength  i  =  1.06  urn  as  a  function 
of  the  intensity  parameter  v  for  different  Q  values. 


It  is  interesting  to  note  that  in  the  relativistic  theory  in 
first  order  of  (hw/mc2)  the  enhancement  R  still  has  the  form 
(7.6).  Just  the  function  n  which  depends  only  on  the  Q  value  of 
the  decay  has  to  be  replaced  by  a  relativistic  expression  [28]. 
This  implies  that  the  field  dependence  of  the  enhancement  R  is 


SUMMARY 


We  have  shown  that  the  modification  of  the  beta-decay  life¬ 
time  due  to  an  external  electromagnetic  field  is  a  pure  quantum 
effect.  For  presently  feasible  laser  fields  the  interaction 
between  the  emitted  electron  and  the  field  is  well  described  by 
the  classical  limit  h  -*•  0,  so  that  the  nuclear  lifetime  is  not 
affected.  In  order  to  obtain  appreciable  quantum  effects  in 
nuclear  beta  de  ays  either  ( i )  the  photon  energy  has  to  be  of  the 
order  of  the  Q  value  of  the  decay  [30]  or  (ii)  the  field  strength 
must  be  comparable  with  the  critical  field  strength.  We  have 
furthermore  seen  in  Fig.  4  that  the  modification  of  the  lifetime 
is  to  a  fair  approximation  a  nonrelativistic  quantum  effect. 
Relativistic  quantum  corrections  are  small  for  all  realistic 
nuclei . 

Crucial  for  the  result  that  the  classical  electron-field 
interaction  does  not  give  rise  to  a  change  in  the  nuclear  life¬ 
time,  is  the  interplay  between  the  pA-term  of  the  interaction 
which  is  responsible  for  the  energy  exchange  between  the  electron 
and  the  field,  and  the  £2-term  which  contributes  to  the  effective 
mass.  For  intense  field  problems  the  P^-term  is  very  important 
and  cannot  be  neglected  as  it  is  often  done  in  quantum  optics. 

The  result  found  here  contradicts  directly  Refs.  10  and  14, 
which  predicted  a  considerable  enhancement  of  free  neutron  and 
nuclear  beta  decay  rates  with  lasers  available  at  present  and 
casts  doubt  on  Refs.  13  and  15  in  which  similar  enhancement  for  y 
decay  and  forbidden  beta  decay  are  obtained.  We  can  state  that 
the  numerical  results  in  all  of  these  four  publications  are  (due 
to  completely  different  reasons)  incorrect.  In  Refs.  10  and  14 
the  "enhancement"  of  the  decay  rate  is  due  to  a  careless  handling 
of  sums  over  Bessel  functions  [25]  as  was  discussed  in  these 
lecture  notes.  In  Ref.  13  the  underlying  physical  situation  is 
not  correctly  modeled  [31]  and  the  calculations  in  Ref.  15  are 
misguided  by  a  wrong  interpretation  of  the  non-interacting 
nuclear  wave  function  (2.13)  and  include  an  algebraical  error, 
after  whose  correction  the  proposed  enhancement  of  forbidden 
decay  rates  disappears  [32].  Obviously,  Nature  does  not  want  her 
more  elementary  constituents  to  be  tampered  with  by  inadequate 
means. 
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kinetical  momentum  p  -  eA.  In  order  to  double  check,  we  notice  from  Eq.  (3) 
that  the  form  factor  satisfies  this  requirement.  The  function  E'(u)  can  be 
rewritten  as 


Actually,  we  have  defined  p  as  the  momentum  of  the  electron  outside  of  the 
field  rather  than  the  canonical  momentum.  Let  us  denote  the  latter  for  the 
moment  by  pc .  We  conclude ^ from  the  fact  that  the  external  field  depends  only 
on  u,  that  the  components  pcT  and  pc0  -  pcz  are  conserved.  In  particular, 
they  do  not  change  when  the  field,  as  a  function  of  u,  is  turned  on  or  off. 
Hence,  inasmuch  as  these  components  are  concerned  and  only  these  do  occur,  we 
can  identify  the  asymptotic  with  the  canonical  momentum. 
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Now,  in  order  to  get  the  total  decay  rate,  the  differential  rate  jii|2  has 
to  be  integrated  over  the  momenta  of  the  decay  products.  If  we  change  the 
integration  variable  from  p  to  p  -  eA(u)  and  commute  the  integration  over  p  with 
that  one  over  u,  the  dependence  of  the  external  field  is  entirely  eliminated. 

Let  us  compare  the  preceding  argument  with  our  earlier  note  regarding  the 
decay  of  a  neutral  free  particle.  The  only  difference  lies  in  the  fact  that  in 
the  latter  case  the  integral  over  space  could  be  explicitly  carried  out  yielding 
a  delta- function  as  the  expression  of  momentum  conservation.  In  the  present 
case,  owing  to  the  presence  of  the  nuclear  wave  functions,  this  was  no  longer 
possible  and  we  were  left  with  the  form  factor  (4).  The  decisive  point  of  the 
argument,  the  dependence  of  the  instantaneous  decay  rate  on  the  kinetical 
momenta,  is  exactly  the  same  in  both  cases. 

The  above  procedure  depends  crucially  on  the  fact  that  we  were  able  to 
write  the  total  rate,  within  the  quasiclassical  approximation,  as  an  integral 
over  an  instantaneous  rate,  as  represented  in  Ea.  (6).  We  expect  that  this 
approximation  will  break  down  as  soon  as  the  external  field  is  so  strong  that 
the  change  in  the  kinetical  momentum  during  the  decay  time  t  is  comparable  with 
the  electron  rest  mass,  viz. 

A(|p  -  eA(t)  | }  -v  e  j  E  (t )  j  It  •v  m  .  (8) 

From  the  uncertainty  relation  we  estimate  It  'v  h/E0  with  Ed  =  -  Efthe  nuclear 

energy  release,  so  that  we  expect  cur  previous  considerations  to  apply  as  long 
as 


In  fact,  calculated  enhancements  of  allowed  [9]  as  well  as  forbidden  [6]  nuc 
beta-decay  are  proportional  to  (m/E0)  3(  |£|/E  ,r^t)2  in  rough  agreement  with  t 
hand-waving  estimate  (9).  To  give  an  explicit  example  we  obtain  for  the  fir 
forbidden  decay  of  90Sr  in  the  presence  of  a  radiofrequency  field  with  \  [;  /•; 
=  6x10  13  and  a  wave  length  of  X  =  100m  an  enhancement  of  2x10* 24 ,  i.e.  zer 
for  all  practical  purposes  in  agreement  with  our  above  argument.  In  ccntras 
in  Ref.  2  a  value  of  the  order  of  unity  is  given. 

To  summarize,  we  have  shown  via  a  very  simple,  general  and  fully 

relativistic  argument  that  nuclear  b.  a-decay,  forbidden  or  not,  cannot  be 

enhanced  by  experimentally  feasible  radiation  fields.  The  situation  is 

particularly  disillusioning  for  the  radio  frequency  fields  with 

|E|/E  ^  10~12,  which  were  considered  in  2tef.  2.  It  is  different  for 

crit 

high- intensity  x-ray  lasers  which  are  nearly  resonant  with  roc  lean  eneray  le 
differences  [10],  since  in  this  case  the  interaction  between  the  nucleus  and 
field  is  of  paramount  importance. 
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